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XI. On the Forces, Stresses, and Fvuxes of Energy in the Electromagnetic Field. 

By Oliver Heaviside, F.R.S. 

Received June 9, — Read June 18, 1891.* 

General Remarks, especially on the Mux of Energy. 

§ 1. The remarkable experimental work of late years has inaugurated a new era in 
the development of the Faraday-Maxwellian theory of the ether, considered as the 
primary medium concerned in electrical phenomena — electric, magnetic, and electro- 
magnetic. Maxwell's theory is no longer entirely a paper theory, bristling with 
unproved possibilities. The reality of electromagnetic waves has been thoroughly 
demonstrated by the experiments of Hertz and Lodge, Fitzgerald and Trouton, 
J. J. Thomson, and others ; and it appears to follow that, although Maxwell's theory 
may not be fully correct, even as regards the ether (as it is certainly not fully 
comprehensive as regards material bodies), yet the true theory must be one of the 
same type, and may probably be merely an extended form of Maxwell's. 

No excuse is therefore now needed for investigations tending to exhibit and 
elucidate this theory, or to extend it, even though they be of a very abstract nature. 
Every part of so important a theory deserves to be thoroughly examined, if only to 
see what is in it, and to take note of its unintelligible parts, with a view to their 
future explanation or elimination. 

§ 2. Perhaps the simplest view to take of the medium which plays such a necessary 
part, as the recipient of energy, in this theory, is to regard it as continuously filling all 
space, and possessing the mobility of a fluid rather than the rigidity of a solid. If 
whatever possess the property of inertia be matter, then the medium is a form of 
matter. But away from ordinary matter it is, for obvious reasons, best to call it as 
usual by a separate. name, the ether. Now, a really difficult and highly speculative 
question, at present, is the connection between matter (in the ordinary sense) and 
ether. When the medium transmitting the electrical disturbances consists of ether 
and matter, do they move together, or does the matter only partially carry forward 
the ether which immediately surrounds it ? Optical reasons may lead us to conclude, 
though only tentatively, that the latter may be the case; but at present, for the 
purpose of fixing the data, and in the pursuit of investigations not having specially 

* Typographical troubles have delayed the publication of this paper. The footnotes are of date 
May 11, 1892. 
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optical bearing, it is convenient to assume that the matter and the ether in contact 
with it move together. This is the working hypothesis made by H. Hertz in his 
recent treatment of the electrodynamics of moving bodies; it is, in fact, what we 
tacitly assume in a straightforward and consistent working out of Maxwell's 
principles without any plainly- expressed statement on the question of the relative 
motion of matter and ether ; for the part played in Maxwell's theory by matter is 
merely (and, of course, roughly) formularised by supposing that it causes the etherial 
constants to take different values, whilst introducing new properties, that of dissipating 
energy being the most prominent and important. We may, therefore, think of merely 
one medium, the most of which is uniform (the ether), whilst certain portions (matter 
as well) have different powers of supporting electric displacement and magnetic induc- 
tion from the rest, as well as a host of additional properties ; and of these we can 
include the power of supporting conduction current with dissipation of energy according 
to Joule's law, the change from isotropy to eolotropy in respect to the distribution 
of the several fluxes, the presence of intrinsic sources of energy, &c. # 

§ 3. We do not in any way form the equations of motion of such a medium, even 
as regards the uniform simple ether, away from gross matter ; we have only to discuss 
it as regards the electric and magnetic fluxes it supports, and the stresses and fluxes 
of energy thereby necessitated. First, we suppose the medium to be stationary, and 
examine the flux of electromagnetic energy. This is the Poyntino flux of energy. 
Next we set the medium into motion of an unrestricted kind. We have now neces- 
sarily a convection of the electric and magnetic energy, as well as the Poynting 
flux. Thirdly, there must be a similar convection of the kinetic energy, &c, of the 
translation al motion ; and fourthly, since the motion of the medium involves the 
working of ordinary (Newtonian) force, there is associated with the previous a flux of 
energy due to the activity of the corresponding stress. The question is therefore a 
complex one, for we have to properly fit together these various fluxes of energy in 
harmony with the electromagnetic equations. A side issue is the determination 
of the proper measure of the activity of intrinsic forces, when the medium moves ; in 
another form, it is the determination of the proper meaning of " true current " in 
Maxwell's sense. 

§ 4. The only general principle that we can bring to our assistance in interpreting 
electromagnetic results relating to activity and flux of energy, is that of the per- 

* Perhaps it is best to say as little as possible at present about the connection between matter and 
ether, but to take the electromagnetic equations in an abstract manner. This will leave us greater 
freedom for f uture modifications without contradiction. There are, also, cases in which it is obviously 
impossible to suppose that matter in bulk carries on with it the ether in bulk which permeates it. 
Either, then, the mathematical machinery must work between the molecules ; or else, we must make such 
alterations in the equations referring to bulk as will be practically equivalent in effect. For example, 
the motional magnetic force VDq of equations (88), (92), (93) may be modified either in q or in D, by 
use of a smaller effective velocity q, or by the substitution in D or cE of a modified reckoning of c for 
the effective permittivity. 
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sistence of energy. But it would be quite inadequate in its older sense referring to 
integral amounts ; the definite localisation by Maxwell, of electric and magnetic 
energy, and of its waste, necessitates the similar localisation of sources of energy ; 
and in the consideration of the supply of energy at certain places, combined with the 
continuous transmission of electrical disturbances, and therefore of the associated 
energy, the idea of a flux of energy through space, and therefore of the continuity of 
energy in space and in time, becomes forced upon us as a simple, useful, and necessary 
principle, which cannot be avoided. 

When energy goes from place to place, it traverses the intermediate space. Only 
by the use of this principle can we safely derive the electromagnetic stress from the 
equations of the field expressing the two laws of circuitation of the electric and 
magnetic forces ; and this again becomes permissible only by the postulation of the 
definite localisation of the electric and magnetic energies. But we need not go so far 
as to assume the objectivity of energy. This is an exceedingly difficult notion, and 
seems to be rendered inadmissible by the mere fact of the relativity of motion, on 
which kinetic energy depends. We cannot, therefore, definitely individualise energy 
in the same way as is done with matter. 

If p be the density of a quantity whose total amount is invariable, and which can 

change its distribution continuously, by actual motion from place to place, its equation 

of continuity is 

conv q/> = />, (1) 

where q is its velocity, and qp the flux of p. That is, the convergence of the flux of 
p equals the rate of increase of its density. Here p may be the density of matter. 
But it does not appear that we can apply the same method of representation to the 
flux of energy. We may, indeed, write 

conv X = T, . . . . (2) 

if X be the flux of energy from all causes, and T the density of localisable energy. 
But the assumption X = Tq would involve the assumption that T moved about like 
matter, with a definite velocity. A part of T may, indeed, do this, viz., when it is 
confined to, and is carried by matter (or ether) ; thus we may write 

conv (qT + X) = T, . . . / (3) 

where T is energy which is simply carried, whilst X is the total flux of energy from 
other sources, and which we cannot symbolise in the form Tq; the energy which 
comes to us from the Sun, for example, or radiated energy. It is, again, often 
impossible to carry out the principle in this form, from a want of knowledge of how 
energy gets to a certain place. This is, for example, particularly evident in the case of 
gravitational energy, the distribution of which, before it is communicated to matter, 
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increasing its kinetic energy, is highly speculative. If it come from the ether (and 
where else can it come from ?), it should be possible to symbolise this in X, if not in 
qT ; but in default of a knowledge of its distribution in the ether, we cannot do so, 
and must therefore turn the equation of continuity into 

S + conv (qT + X) = T, ............' (4) 

where S indicates the rate of supply of energy per unit volume from the gravitational 
source, whatever that may be. A similar form is convenient in the case of intrinsic 
stores of energy, which we have reason to believe are positioned within the element 
of volume concerned, as when heat gives rise to thermoelectric force. Then S is the 
activity of the intrinsic sources. Then again, in special applications, T is conveniently 
divisible into different kinds of energy, potential and kinetic. Energy which is 
dissipated or wasted comes under the same category, because it may either be regarded 
as stored, though irrecoverably, or passed out of existence, so far as any immediate 
useful purpose is performed. Thus we have as a standard practical form of the 
equation of continuity of energy referred to the unit volume, 

S + conv {X + q (U + T)} = Q + TJ + T, . . . . . . . . . (5) 

where S is the energy supply from intrinsic sources, U potential energy and T kinetic 
energy of localisable kinds, q (TJ + T) its convective flux, Q the rate of waste of 
energy, and X the flux of energy other than convective, e.g., that due to stresses in 
the medium and representing their activity. In the electromagnetic application we 
shall see that TJ and T must split into two kinds, and so must X, because there is a 
flux of energy even when the medium is at rest. 

§ 5. Sometimes we meet with cases in which the flux of energy is either wholly or 
partly of a circuital character. There is nothing essentially peculiar to electromagnetic 
problems in this strange and apparently useless result. The electromagnetic instances 
are paralleled by similar instances in ordinary mechanical science, when a body is 
in motion and is also strained, especially if it be in rotation. This result is a necessary 
consequence of our ways of reckoning the activity of forces and of stresses, and serves 
to still further cast doubt upon the " thinginess " of energy. At the same time, the 
flux of energy is going on all around us, just as certainly as the flux of matter, and 
it is impossible to avoid the idea ; we should, therefore, make use of it and formularise 
it whenever and as long as it is found to be useful, in spite of the occasional failure to 
obtain readily understandable results. 

The idea of the flux of energy, apart from the conservation of energy, is by no 
means a new one. Had gravitational energy been less obscure than it is, it might 
have found explicit statement long ago. Professor Poynting* brought the principle 

* Potnting, 'Phil Trans.,' 1884. 
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into prominence in 1884, by making use of it to determine the electromagnetic flux of 
energy. Professor Lodge # gave very distinct and emphatic expression of the principle 
generally, apart from its electromagnetic aspect, in 1885, and pointed out how much 
more simple and satisfactory it makes the principle of the conservation of energy 
become. So it would, indeed, could we only understand gravitational energy ; but in 
that, and similar respects, it is a matter of faith only. But Professor Lodge attached, 
I think, too much importance to the identity of energy, as well as to another principle 
he enunciated, that energy cannot be transferred without being transformed, and 
conversely ; the transformation being from potential to kinetic energy or conversely. 
This obviously cannot apply to the convection of energy, which is a true flux of 
energy ; nor does it seem to apply to cases of wave motion in which the energy, 
potential and kinetic, of the disturbance, is transferred through a medium unchanged 
in relative distribution, simply because the disturbance itself travels without change 
of type ; though, it may be that in the unexpressed internal actions associated with 
the wave propagation there might be found a better application. 

It is impossible that the ether can be fully represented, even merely in its trans- 
missive functions, by the electromagnetic equations. Gravity is left out in the cold ; 
and although it is convenient to ignore this fact, it may be sometimes usefully 
remembered, even in special electromagnetic work ; for, if a medium have to contain 
and transmit gravitational energy as well as electromagnetic, the proper system of 
equations should show this, and, therefore, include the electromagnetic. It seems, 
therefore, not unlikely that in discussing purely electromagnetic speculations, one 
may be within a stone's throw of the explanation of gravitation all the time. The 
consummation would be a really substantial advance in scientific knowledge. 

On the Algebra and Analysis of Vectors without Quaternions. Outline of Author's 

System. 

§ 6. The proper language of vectors is the algebra of vectors. It is, therefore, 
quite certain that an extensive use of vector-analysis in mathematical physics 
generally, and in electromagnetism, which is swarming with vectors, in particular, is 
coming and may be near at hand. It has, in my opinion, been retarded by the want of 
special, treatises on vector analysis adapted for use in mathematical physics, Professor 
Tait's well-known profound treatise being, as its name indicates, a treatise on 
Quaternions. I have not found the Hamilton-Tait notation of vector operations 
convenient, and have employed, for some years past, a simpler system. It is not, 
however, entirely a question of notation that is concerned. I reject the quaternionic 
basis of vector-analysis. The anti-quaternionic argument has been recently ably 
stated by Professor Willard Gibbs.1" He distinctly separates this from the question 

* Lodge, ' Phil. Mag.,' June, 1885, " On the Identity of Energy." 

t Professor Gibbs's letters will be found in 'Nature,' vol. 43, p. 511, and vol. 44, p. 79 ; and Professor 
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of notation, and this may be considered fortunate, for whilst I can fully appreciate and 
(from practical experience) endorse the anti-quaternionic argument, I am unable to 
appreciate his notation, and think that of Hamilton and Tait is, in some respects, 
preferable, though very inconvenient in others. 

In Hamilton's system the quaternion is the fundamental idea, and everything 
revolves round it. This is exceedingly unfortunate, as it renders the establishment 
of the algebra of vectors without metaphysics a very difficult matter, and in its 
application to mathematical analysis there is a tendency for the algebra to get more 
and more complex as the ideas concerned get simpler, and the quaternionic basis forms 
a real difficulty of a substantial kind in attempting to work, in harmony with ordinary 
Cartesian methods. 

Now, I can confidently recommend, as a really practical working system, the 
modification I have made. It has many advantages, and not the least amongst them is 
the fact that the quaternion does not appear in it at all (though it may, without 
much advantage, be brought in sometimes), and also that the notation is arranged so 
as to harmonise with Cartesian mathematics. It rests entirely upon a few definitions, 
and may be regarded (from one point of view) as a systematically abbreviated 
Cartesian method of investigation, and be understood and practically used by any one 
accustomed to Cartesians, without any study of the difficult science of Quaternions. 
It is simply the elements of Quaternions without the quaternions, with the notation 
simplified to the uttermost, and with the very inconvenient minus sign before scalar 
products done aw 7 ay with/" 

Tait's in vol. 43, pp. 535, 608. This rather one-sided discussion, arose out of Professor Tait stigmatising 
Professor GIibbs as " a retarder of quaternionic progress." This may be very true ; but Professor GibbS 
is anything but a retarder of progress in vector analysis and its application to physics. 

# §§ 7, 8, 9 contain an introduction to vector-analysis (without the quaternion), which is sufficient 
for the purposes of the present paper, and, I may add, for general use in mathematical physics. It is an 
expansion of that given in my paper " On the Electromagnetic Wave Surface," ' Phil. Mag.,' June, 1885. 
The algebra and notation are substantially those employed in all my papers, especially in " Electromag- 
netic Induction and its Propagation, " ' The Electrician,' 1885. 

Professor Gibbs's vectorial work is scarcely known, and deserves to be well known. In June, 1888, I 
received from him a little book of 85 pages, bearing the singular imprint Not Published. Newhaven, 
1881-4. It is indeed odd that the author should not have published what he had been at the trouble of 
having printed. His treatment of the linear vector operator is specially deserving of notice. Although 
" for the use of students in physics," I am bound to say that I think the work much too condensed for a 
first introduction to the subject. 

In 'The Electrician' for Nov. 13, 1891, p. 27, I commenced a few articles on elementary vector- 
algebra and analysis, specially meant to explain to readers of my papers how to work vectors. I am 
given to understand that the earlier ones, on the algebra, were much appreciated; the later ones, 
however, are found difficult. But the vector-algebra is identically the same in both, and is of quite a 
rudimentary kind. The difference is, that the later ones are concerned with analysis, with varying 
vectors; it is the same as the difference between common algebra and differential calculus. The 
difficulty, whether real or not, does not indicate any difficulty in the vector-algebra. I mention this on 
account of the great prejudice which exists against vector-algebra. 
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§ 7. Quantities being divided into scalars and vectors, I denote the scalars, as 
usual, by ordinary letters, and put the vectors in the plain black type, known, I 
believe, as Clarendon type, rejecting Maxwell's German letters on account of their 
being hard to read. A special type is certainly not essential, but it facilitates the 
reading of printed complex vector investigations to be able to see at a glance which 
quantities are scalars and which are vectors, and eases the strain on the .memory. 
But in MS. work there is no occasion for specially formed letters. 

Thus A stands for a vector. The tensor of a vector may be denoted by the same 
letter plain ; thus A is the tensor of A. (In MS. the tensor is A .) Its rectangular 
scalar components are A 1? A 2? A 3 . A unit vector parallel to A may be denoted by 
A l5 so that A = AA X . But little things of this sort are very much matters of taste. 
What is important is to avoid as far as possible the use of letter prefixes, which, 
when they come two (or even three) together, as in Quaternions, are very confusing. 

The scalar product of a pair of vectors A and B is denoted by AB, and is defined 
to be 

A 

AB = A L B X + A 3 B 3 + A 3 B 3 = AB cos AB = BA . (6) 

The addition of vectors being as in the polygon of displacements, or velocities, or 
forces ; i.e., such that the vector length of any closed circuit is zero ; either of the 
vectors A and B may be split into the sum of any number of others, and the multi- 
plication of the two sums to form AB is done as in common algebra ; thus 

(a + b) (c + d) = ac -f ad -f be + bd, 1 

— ca -f da 4- cb -f- db. J v 

If N be a unit vector, NN or K 2 = 1 ; similarly A 2 = A 2 for any vector. 

The reciprocal of a vector A has the same direction ; its tensor is the reciprocal of 
the tensor of A. Thus 

A 

AAA 
AB 1 = B" 1 A = ~ = ~ cos AB (8) 



and 



The vector product of a pair of vectors is denoted by VAB, and is defined to be 

A 

the vector whose tensor is AB sin AB, and whose direction is perpendicular to 
the plane of A and B. Or 

VAB = i (A 2 B 8 - A 3 B 3 ) + j (A 3 B X - A X B 3 ) + k (A X B 3 - A 2 B X ) - - YBA, ... (9) 

where i, j, k, are any three mutually rectangular unit vectors. The tensor of VAB is 
Y AB ; or 

V AB = AB sin AB. . . (10) 
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Its components are iVAB, jVAB, kVAB. 

In accordance with the definitions of the scalar and vector products, we have 



j* = l, 

ik = 0, 



k 2 ^ 1 
ki=0 



i, Vki = j 



> 



(11) 



i 2 = l, 

ij = 0, 
Vij = k, 

and from these we prove at once that 

Y (a + b) (c + d) = Vac + Yad + Vbc + Ybd, 

and so on, for any number of component vectors. The order of the letters in each 
product has to be preserved, since Vab = — Vba, 
Two very useful formulae of transformation are 



and 



= BVCA = GYAE 
A, (B 2 C 3 - B 3 C 2 ) + A 3 (BA - BA) + As (*A - B^) ; 

YAYBC = B.CA - CAB, 1 
= B(CA) ~C(AB). 



» a 



» * 



# a * 



(12) 



. . (13) 



Here the dots, or the brackets in the alternative notation, merely act as separators, 
separating the scalar products CA and AB from the vectors they multiply. A space 
would be equivalent, but would be obviously unpractical. 



A . 



A 



As — is a scalar product, so in harmony therewith, there is the vector product V - . 
Since VAB = — VBA, it is now necessary to make a convention as to whether the 

A 

denominator comes first or last in V-. Say therefore, VAB" 1 , Its tensor is 



B 



v A _ A . A 
Vq^d — o sin AB, 



a a s 



• (14) 



§ 8. Differentiation of vectors, and of scalar and vector functions of vectors with 
respect to scalar variables is done as usual. Thus, 



d 

dt 



a a i) 9 

A = iA T -f j A 3 + kA 3 . 



AB = AB + BA. 



•^ 



>■ • 



• * 



n AYBC = AYBC + AYBC -h AYBC. 



at 



. . (15) 



j 



The same applies with complex scalar differentiators, e.g., with the differentiator 



d ___ d 

dt~~~ dt 



4-qV, 
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used when a moving particle is followed, q being its velocity. Thus, 

3 9b 3a 

dt dt dt 

= AB + BA + A.qV.B + B.qVJL ......... (16) 

Here qV is a scalar differentiator given by 

q V = ft ^ + ffa ^ + 2 S £;, .......... (17) 

so that A.qV.B is the scalar product of A and the vector qV.B ; the dots here again 
act essentially as separators. Otherwise, we may write it A (qV) B. 
The fictitious vector V given by 

V = iV 1+ jV 3 + kV 3 = i|;+j| + k| (18) 

is very important. Physical mathematics is very largely the mathematics of V. The 
name Nabla seems, therefore, ludicrously inefficient. In virtue of i, j, k, the operator 
V behaves as a vector. It also, of course, differentiates what follows it. 
Acting on a scalar P, the result is the vector 

VP = iV x P + jVoP + kV 3 P : , (19) 

the vector rate of increase of P with length. 

If it act on a vector A, there is first the scalar product 

VA = V 1 A 1 + V 3 A 3 + V 3 A 3 = div A, (20) 

or the divergence of A. Regarding a vector as a flux, the divergence of a vector is 
the amount leaving the unit volume. 
The vector product VVA is 

YVA - i (V 3 A 3 - V 8 A 3 ) + j (VA - VAs) + k (V x A 3 - V 3 A X ) 5 

== curl A (21) 

The line-integral of A round a unit area equals the component of the curl of A 
perpendicular to the area. 

We may also have the scalar and vector products NV and VNV, where the vector 
N is not differentiated. These operators, of course, require a function to follow them 
on which to operate ; the previous qV.A of (16) illustrates. 

The Laplacean operator is the scalar product V 2 or VV ; or 

V 3 = Vi 2 + Vjj 3 + V 3 * ; . (22) 
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and an example of (13) is 

VVVVA = V.VA - V 2 A 3 
or 

cuvPA = V div A — V 3 A, ........... (23) 

which is an important formula. 

Other important formulae are the next three, 

div PA = P div A + AV.P, ........... (24) 

P being scalar. Here note tha.t AV.P and AVP (the latter being the scalar product 
of A and VP) are identical. This is not true when for P we substitute a vector. 
Also 

div YAB = B curl A •— A curl B ; . . . . (25) 

which is an. example of (12), noting that both A and B have to be differentiated, 

jcx.ncl 

curl YAB = BV.A + A div B — AVJ3 — B div A, ....... (26) 

This is an example of (13). 

§ 9o When one vector D is a linear function of another vector E, that is, con- 
nected by equations of the form 

Di = % E i + %E 3 + c 13 E 3 r 

D 3 = c 21 E 1 + c 33 E 3 + c 23 E 3 , V, ........... (27) 

D 8 = c^ + c 82 E 2 + c 33 E 3 J 

in terms of the rectangular components, we denote this simply by 
where c is the linear operator. The conjugate function is given by 



it j 



where D' is got from D by exchanging c 12 and c 31 , &c« Should the nine coefficients 
reduce to six by c 13 = c 2 i> & c «> D and D' are identical, or D is a self-conjugate or sym- 
metrical linear function of B. 

But, in general, it is the sum of D and I)' which is a symmetrical function of E, 
and the difference is a simple vector-product. Thus 

c E + V€E, L 30 



* « « 






where c is a self-conjugate operator, and € is the vector given by 
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2 2 2 ........ . 

The important characteristic of a self-conjugate operator is 

where E x and E 3 are any two E's, and D l9 D 2 the corresponding D's. But when 
there is not symmetry, the corresponding property is 

Ei Da = E 9 D -i, 

or E lC E 3 = E 3 c'E r r v J 

Of these operators we have three or four in electromagnetism connecting forces and 
fluxes, and three more connected with the stresses and strains concerned. As it 
seems impossible to avoid the consideration of rotational stresses in electromagnetism, 
and these are not usually considered in works on elasticity, it will be desirable to 
briefly note their peculiarities here, rather than later on. 



On Stresses, irrotational and rotational, and their Activities. 

§ 10* Let P N be the vector stress on the N plane, or the plane whose unit normal 
is N. It is a linear function of N. This will fully specify the stress on any plane. 
Thus, if P 1? P 3 , P 3 are the stresses on the i, j, k planes, we shall have 

*i = iPn + jP 19 + kP ]8 r 

*3 = 1 ^ ) 2i + JP22 + kP & > ........... (34) 

P 3 = iP 31 + J'P 3 2 + kPffl- J 

Let, also, Q N be the conjugate stress ; then, similarly, 

i = iP n +jPgi + kPgi,- 

s = iP 13 + jP 23 + kP 33 , V . . . . ,. . . , , , , (35) 



$3 = iPl3 + J* P 23 + k P 3 3- 



J 



Half the sum of the stresses JP N and Q N is an ordinary irrotational stress ; so that 

3P N = $JK + VgN, 1 

a N = o K-v€N,J (36) 

where (f> is self-conjugate, and 

MDCCCXCIT.— A. 3 k 
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2€ = i (P 23 — P 33 ) + j (P 31 — P 13 ) + k (Pis — Psi)» ....... (37) 



Here 2€ is the torque per unit volume arising from the stress 

The translations] force, F, per unit volume is (by inspection of a unit cube) 

p = V^ + V 3 P 3 + V 3 P 3 ■.....'........ (38) 

= i div ftj + j div Q 3 -f k div €t 3 ;......... . (39) 

or, in terms of the self-conjugate stress and the torque, 

F = (i div o i -f j div o j + k div o k) — curl € ? . . * . . . . (40) 

where — curl € is the translational force due to the rotational stress alone, as in Sir 
W. Thomson's latest theory of the mechanics of an " ether e " # 

Next, let N be the unit normal drawn outward from any closed surface. Then 



N 



CtOQB»1»aO«98« 



(41) 



where the left summation extends over the surface and the right summation through- 
out the enclosed region. For 

P N = N 1 P 1 + N a P 2 + N 3 P 8 

=:i.E , a l -f-j.E , a 3 + k.m 3 ; ........... (42) 

so the well-known theorem of divergence gives immediately, by (39), 

2P N = 2 (i div ft x + j div a 3 + k div a 3 ) = 2 P. ....... (43) 

Next, as regards the equivalence of rotational effect of the surface-stress to that of 
the internal forces and torques. Let r be the vector distance from any fixed origin. 
Then VrF is the vector moment of a force^ F, at the end of the arm r. Another 
(not so immediate) application of the divergence theorem gives 

2* VrPjj =s & Vr.fi ~j- ^2€? « • • • * ° • * s ° {jb^j 

Thus, any distribution of stress, whether rotational or irrotational, may be regarded as 
in equilibrium. Given any stress in a body, terminating at its boundary, the body 
will be in equilibrium both as regards translation and rotation. Of course, the 
boundary discontinuity in the stress has to be reckoned as the equivalent of internal 
divergence in the appropriate manner. Or, more simply, let the stress fall off 
continuously from the finite internal stress to zero through a thin surface-layer. We 

* 'Mathematical and Physical Papers/ vol. 3, Art. 99, p. 436. 
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then have a distribution of forces and torques in the surface- layer which equilibrate 
the internal forces and torques. 

To illustrate; we know that Maxwell arrived at a peculiar stress, compounded of 
a tension parallel to a certain direction, and an equal lateral pressure, which would 
account for the mechanical actions apparent between electrified bodies, and endea- 
voured similarly to determine the stress in the interior of a magnetised body to 
harmonise with the similar external magnetic stress of the simple type mentioned. 
This stress in a magnetised body I believe to be thoroughly erroneous ; nevertheless, 
so far as accounting for the forcive on a magnetised body is concerned, it will, when 
properly carried out with due attention to surface-discontinuity, answ r er perfectly 
well, not because it is the stress, but because any stress would do the same, the only 
essential feature concerned being the external stress in the air. 

Here we may also note the very powerful nature of the stress-function, considered 
merely as a mathematical engine, apart from physical reality. For example, we may 
account for the forcive on a magnet in many ways, of which the two most prominent 
are by means of forces on imaginary magnetic matter, and by forces on imaginary 
electric currents, in the magnet and on its surface. To prove the equivalence of 
these two methods (and the many others) involves very complex surface- and volume- 
integrations and transformations in the general case, w T hich may be all avoided by the 
use of the stress-function instead of the forces. 

§ 11. Next as regards the activity of the stress P N and the equivalent translational, 
distortional, and rotational activities. The activity of P N is P N q per unit area, if q 
be the velocity. Here 

P»q = ft-JSrO! + ga-NQa + ^BTOb, (45) 

by (42) ; or, re-arranging, 

P N q = H (ftttx + q%% + g 3 a 3 ) = N20g, 

where Q^ is the conjugate stress on the q plane,. That is, qQ, g or 2Qg is the 
negative of the vector flux of energy expressing the stress-activity. For we choose 
P NN so as to mean a pull when it is positive, and when the stress P N works in the 
same sense with q energy is transferred against the motion, to the matter which is 
pulled. 

The convergence of the energy-flux, or the divergence of qQ, q , is therefore the 
activity per unit volume. Thus 

div (0^ + O^ + a 3 2 3 ) = q (i div ^ + j div % + k div Qg) + (O^i + GgVft + a s Vg 3 ) . (47) 

- q (VA + V 3 P 3 + V 3 P 3 ) + P x V x q + P 3 V 3 q + P 3 V 3 q .... (48) 

where the first form (47) is generally most useful. Or 
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div 2 Q,q = Fq + 2 GVg ; ........... (49) 

where the first term on the right is the translation^] activity, and the rest is the sum 
of the distortional and rotational activities. To separate the latter introduce the 
strain velocity vectors (analogous to P 1? P 2 , P 3 ) 

Pi = i(Vff 1 + V 1 4), P 2 = ir( v 22 + V 8 q), p 3 = | (Vg 3 + V 3 q) ; . . . . (50) 

and generally 

p N = i (V.qN + NV.q). ............ (51) 

Using these we obtain 

= 2 ftp + | OjVi curl <1 + | ft 3 Vj curl q + £ ft 3 Vk curl <i 

«*— — ^m/ HCy JLr J " %j# L/ Li. X X v| «e e 9 g <s e o 9 « e . * e » » e q « e o \kjjmlj 

Thus 2Qp is the distortional activity and € curl q the rotational activity. But 
since the distortion and the rotation are quite independent, we may put 2Pp for the 
distortional activity ; or else use the self-conjugate stress, and wxite it ^ 2 (P + Q) p. 

§ 12. In an ordinary " elastic solid," when isotropic, there is elastic resistance to 
compression and to distortion. We may also imaginably have elastic resistance to 
translation and to rotation; nor is there, so far as the mathematics is concerned, any 
reason for excluding dissipative resistance to translation, distortion, and rotation ; 
and kinetic energy may be associated with all three as well, instead of with the 
translation alone, as in the ordinary elastic solid. 

Considering only three elastic moduli, we have the old h and n of Thomson and 
Tait (resistance to compression and rigidity), and a new coefficient, say n v such that 

^J _ ™ %-i CUrj, Aj-j 8(j«i9«o»c«e«Oe \^OtjJ 

if D be the displacement and 2€ the torque^ as before 
The stress on the i plane (any plane) is 

P 1 = w (VDj + V X D) + i (ft - | n) div D + ^ V curl D.i 

=5= (n + %) V{D + (n — n x ) VJ) 1 + (ft — f ») i div D ; . . . . . . (54) 

and its conjugate is 

a x = » (VDjl + V X D) + i (ft - f n) div 3D - % ( V X D - VD^ 

= (% — w x ) V^ -h (w ■+■ ^i) VD X -f i (ft — f ft) div D j ».»».. (55) 
from which 

\ =s div C^ s= (w — ^ + ft — f ») ^i div D + (*& + %) V 2 D X o * . . , (56) 



is the i component of the translational force ; the complete force F is therefore 
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F = O + n x ) V Z B + (k -f %n - nj V div D ; . . . . . / . . (57) 

or, in another form, if 

P = - A; div B, 

P being the isotropic pressure, 

F = - VP + n (V 3 B + J V div B) - n Y curl 3 B, (58) 

remembering (23) and (53). 

We see that in (57) the term involving div D may vanish in a compressible solid 
by the relation n x = k + | ft ; this makes 

^ -f- ^ = ft + |- n, n x — w = k — %n, . . . . . . , . . (59) 

which are the moduli, longitudinal and lateral, of a simple longitudinal strain ; that 
is, multiplied by the extension, they give the longitudinal traction, and the lateral 
traction required to prevent lateral contraction. 

The activity per unit volume, other than translational, is 

2aV 2 = (w - n L ) (^D.Vft + V 3 B.Vg 3 + V 3 B.Vg 3 ) 
+■ (n + n{) (VB i ,Wq l + VD 2 .Vg 3 + VD 3 .V? 3 ) 

= n (V^.Vg, + V 2 D.V ?2 + V 8 D.Vg s + V^V^ + VD 2 V? 2 + VD 3 V 2s ) 

4- (A; — f w) div D div qfwj curl B curl q ;.,,.. . ..... . (60) 

or, which is the same, 

2QV 2 = ~ [\k (div B) 3 + fa (curl B) 3 

+ M ( VD i) 2 + (VD 3 ) 2 + (VD 3 )2 + VD^D + VD 3 .V 3 B + VD 3 .V 3 B j - ^ (div D)*], . (6i) 

where the quantity in square brackets is the potential energy of an infinitesimal 
distortion and rotation. The italicised reservation appears to be necessary, as we 
shall see from the equation of activity later, that the convection of the potential 
energy destroys the completeness of the statement 

2qv 2 == TJ, 
if U be the potential energy. 

In an elastic solid of the ordinary kind, with n x =s 0, we have 

f N = n (2 curl YBU + VSf curl D), 
1? == — n curl 3 Bo 



In the case of a medium in which n is £ero but n ± finite (Sir W, Thomson's 
rotational ether) 5 
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P N = %! V curl D.N, "1 

F = - ^ curl* D. J ^ 

Thirdly, if we have both Jc = — f n and n = b 1? then 

P H = 2n curl VDN, 

o wt* >(€=*»ourlD), ...... (64) 

F = — 2w curl 3 3D, 

i.e., the sums of the previous two stresses and forces, 

§13. As already observed^ the vector flux of energy, due to the stress, is 

s 

— 2% =s — Q^q = — (tttfi + O3& + % 3 )« *.......* (65) 

Besides this, there is the flux of energy 

4 (U + T) 

by convection^ where U is potential and T kinetic energy, Therefore, 

W = 4(U + T) -~20g . . . . * . * . . ", . . . (66) 

represents the complete energy flux, so far as the stress and motion are concerned. 
Its convergence increases the potential energy, the kinetic energy, or is dissipated. 
But if there be an impressed translational force f its activity is % This supply of 
energy is independent of the convergence of W. Hence 

f4 + Q + U==T + div[q(U + T)-20g] ........ (67) 

is the equation of activity. 

But this splits into two parts at least, For (67) is the same as 

(f+F)cL + 2aVg = Q + U + T + div4(tJ + T), ,..,... (68) 
and the translational portion may be removed altogether. That is, 

(f+F)<l=Q + U + T + dIvct(U + T ) ? a ....... (69) 

if the quantities with the zero suffix are only translationally involved. For 
example, if 

ut 

as in fluid motion, without frictional or elastic forces associated with the translation, 
then 
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(f + F) q = P (L ~ = T + div qT, . (71) 

if T = ^pq 2 , the kinetic energy per unit volume. The complete form (69) comes 
in by the addition of elastic and frictional resisting forces. So deducting (69) from 
(68) there is left 

2<lVg = Qi + Ui + Tj + div q (U L + T x ), ......... (72) 

where the quantities with suffix unity are connected with the distortion and the 
rotation, and there may plainly be two sets of dissipative terms, and of energy (stored) 
terms. Thus the relation 



€ 



= ( w 3 +%~ H- w 3 — ) curl D .......... (73) 



will bring in dissipation and kinetic energy, as well as the former potential energy of 
rotation associated with n v 

That there can be dissipative terms associated with the distortion is also clear 
enough, remembering Stokes's theory of a viscous fluid. Thus, for simplicity, do away 
with the rotating stress, by putting € = 0, making P N and Q N identical. Then take 
the stress on the i plane to be given by 

p ^( w+ ^l + "S) (VD ^ v ^- i { p+ *( w + ^, + "S) divD }' • • <**> 

and similarly for any other plane ; where P = — k div D. 

When jx = 0, v = 0, we have the elastic solid with rigidity and compressibility. 
When n = 0, v = 0, we have the viscous fluid of Stokes. When v = only, we 
have a viscous elastic solid, the viscous resistance being purely distortional, and 
proportional to the speed of distortion. But with n, /x, v, all finite, we still further 
associate kinetic energy with the potential energy and dissipation introduced by 
n and /jl. 

We have 

2PVg = Q 8 + U 8 +T 8 

for infinitesimal strains, omitting the effect of convection of energy ; where 

Ta-M-ICdiv^ + VftCVfc^ .... (75) 

= p. [ - | (div q)» + Vft (Vft + V l4 ) + Vg 3 (V? 2 + V 8 q) + Vg 3 (Vg 3 + V 8 q)], .... (76) 



U 3 = |n 
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Observe that T 3 and Q 3 only differ In the exchange of fi to \v ; but U 2 , the potential 
energy, is not the same function of n and D that T 3 is of v and q e But if we take 
h = 0, we produce similarity. An elastic solid having no resistance to compression is 
also one of Sir W. Thomson's ethers. 

When n = 0, /x = 0, v = 5 we come down to the friction less fluid, in which 

f __ yp — n ^L ( 78^ 

dt 

and 

2PVg = — P divq, . (79) 

with the equation of activity 

fq = U + T + div (IT + T + P) q, . . . . . . . . . . (80) 

the only parts of which are not always easy to interpret are the Pq term, and the 
proper measure of U. By analogy, and conformably with more general cases, we 
should take 

p = _ fc div D, and U = i h (div D)\ 

reckoning the expansion or compression from some mean condition. 



The Electromagnetic Equations in a Moving Medium, 

§ 14* The study of the forms of the equation of activity in purely mechanical cases 3 
and the interpretation of the same is useful, because in the electromagnetic problem 
of a moving medium we have still greater generality, and difficulty of safe and sure 
interpretation. To bring it as near to abstract dynamics as possible, all we need say 
regarding the two fluxes, electric displacement B and magnetic induction B, is that 
they are linear functions of the electric force E and magnetic force H, say 



= fii 






where c and /x are linear operators of the symmetrical kind, and that associated with 
them are the stored energies U and T, electric and magnetic respectively (per unit 
volume), given by 

TT — 1 TFTI HP — 1 TTP CR&\ 

In isotropic media c is the permittivity, //, the inductivity* It is unnecessary to 
say more regarding the well-known variability of p and hysteresis than that a magnet 
is here an ideal magnet of constant inductivity* 

As there may be impressed forces, E is divisible into the force of the field and an 
impressed part ; for distinctness, then, the complete E may be called the i6 force of 
the flux ?? De Similarly as regards H and B. 
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There is also waste of energy (in conductors, namely) at the rates 

ic^i ~ iiC, t^2 ~~ UK,., . (8*^) 

where the fluxes C and K are also linear functions of E and H respectively ; thus 

C = frE, K = g1L 9 . \ ......... . (84) 

where, when the force is parallel to the flux, and k is scalar, it is the electric conduc- 
tivity. Its magnetic analogue is g, the magnetic conductivity. That is, a magnetic 
conductor is a (fictitious) body which cannot support magnetic force without con- 
tinuously dissipating energy. 

Electrification is the divergence of the displacement, and its analogue, magnetifica- 
tion, is the divergence of the induction ; thus 

p = div D, <j = div B (85) 

are their volume densities. The quantity a is probably quite fictitious, like K. 

According to Maxwell's doctrine, the true electric current is always circuital, and 
is the sum of the conduction current and the current of displacement, which is the 
time rate of increase of the displacement. But, to preserve circuitality, we must add 
the convection current when electrification is moving, so that the true current 
becomes 

J = C + D + q/>, (86) 

where q is the velocity of the electrification p. Similarly 

G=K + B + q* . (87) 

should be the corresponding magnetic current. 

§ 15, Maxwell's equation of electric current in terms of magnetic force in a 
medium at rest, say, 

curl H^ = C + D, 

where E^ is the force of the field, should be made, using H instead, 

curl (H — h ) = C + D + q/>, 

and here h will be the force of intrinsic magnetisation, such that /xh is the intensity 
of intrinsic magnetisation. But I have shown that when there is motion, another 
impressed term is required, viz., the motional magnetic force 

h = VDq, . . „ (88) 

MDCCCXCII. — A. 3 L 
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making the first circuital law become 

curl (H -\ ~ h) = J = C + 3D + qp. . . , . . . . . . (89) 

Maxwell's other connection to form the equations of propagation is made through 
his vector-potential A and scalar potential V. Finding this method not practically 
workable, and also not sufficiently general, I have introduced instead a companion 
equation to (89) in the form 

— curl (E — e — e) = G = K + B + q*, (90) 

where e expresses intrinsic force, and e is the motional electric force given by 

e = VqB ? ....,....,..., (91) 

which is one of the terms in Maxwell's equation of electromotive force. As for e 0> 
it includes not merely the force of intrinsic electrification, the analogue of intrinsic 
magnetisation, but also the sources of energy, voltaic force, thermoelectric force, &c 

(89) and (90) are thus the working equations, with (88) and (91) in case the 
medium moves ; along with the linear relations before mentioned, and the definitions 
of energy and waste of energy per unit volume. The fictitious K and cr are useful in 
symmetrizing the equations, if for no other purpose. 

Another way of writing the two equations of curl is by removing the e and h 
terms to the right side. Let 

curl h== j, J + j = J G , , 

— curl e = g, G- -f g = " 

Then (89) and (90) may be written 

curl (H — • h ) = J = C + I) -f qp -f j, 
— curl (E — e ) = G = K + B + q* -f g. J 

So far as circuitality of the current goes, the change is needless, and still further 
complicates the make-up of the true current, supposed now to be J . On the other 
hand, it is a simplification on the left side, deriving the current from the force of the 
flux or of the field more simply, 

A question to be settled is whether J or J should be the true current. There 
seems only one crucial test, viz., to find whether e J or e J is the rate of supply of 
energy to the electromagnetic system by an intrinsic force e . This requires, however, 
a full and rigorous examination of all the fluxes of energy concerned. 
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The Electromagnetic Flux of Energy in a stationary Medium. 



§ 16. First let the medium be at rest, giving us the equations 



curl (H - ho) = J = C + D, (94) 

- curl (E - e ) = = H + B (95) 

Multiply (94) by (E — e ), and (95) by (H — h ), and add the results. Thus, 
(E - e () ) J + (H - ho) G = (E - e ) curl (H - h ) - (H - h ) corl (E - e ), 

which, by the formula (25), becomes 

e J + h G = EJ + HG +.div V (E - e ) (H - h ) ; 

or, by the use of (82), (83), 

e J + h G=Q + U + t + clivW, .......... (96) 

where the new vector W is given by 

W = V(E-e )(H-h ). ........... (97) 

The form of (96) is quite explicit, and the interpretation sufficiently clear. The left 

side indicates the rate of supply of energy from intrinsic sources. These (Q + XJ + T) 
shows the rate of waste and of storage of energy in this unit volume. The remainder, 
therefore, indicates the rate at which energy is passed out from the unit volume ; and 
the flux W represents the flux of energy necessitated by the postulated localisation of 
energy and its waste, when E and H are connected in the manner shown by (94) 
and (95). 

There might also be an independent circuital flux of energy, but, being useless, 
it would be superfluous to bring it in. 

The very important formula (97) was first discovered and interpreted by Professor 
Poynting, and independently discovered and interpreted a little later by myself in an 
extended form. It will be observed that in my mode of proof above there is no 
limitation as to homogeneity or isotropy as regards the permittivity, inductivity, and 
conductivity. But c and /x should be symmetrical. On the other hand, k and g do 
not require this limitation in deducing (97).* 

* The method of treating Maxwell's electromagnetic scheme employed in the text (first introduced 
in "Electromagnetic Induction and its Propagation, " < The Electrician,' January 3, 1885, and later) 

3 L 2 
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It is important to recognize that this flux of energy is not dependent upon the 
transitional motion of the medium, for it is assumed explicitly to be at rest. The 
vector W cannot, therefore, be a flux of the kind Q^ before discussed, unless possibly 
it be merely a rotating stress that is concerned. 

The only dynamical analogy with which I am acquainted which seems at all 
satisfactory is that furnished by Sir W. Thomson's theory of a rotational ether. 
Take the case of e = 0, h =: 0, k — Q, g =z 0, and c and fi constants, that is, pure 
ether uncontaminated by ordinary matter. Then 

curl H = eE, ............. (98) 

— curl £i ~ jLtxt, ............. (yyj 

Now, let H be velocity, jjl density; then, by (99), — curl E is the translational 
force due to the stress, which is, therefore, a rotating stress ; thus, 

P N = VE¥, a N = VNE; .......... (100) 

and 2E is the torque. The coefficient c represents the compliancy or reciprocal of 
the quasi-rigidity. The kinetic energy ^/xH 2 represents the magnetic energy, and the 
potential energy of the rotation represents the electric energy; whilst the flux of 

energy is VEIL For the activity of the torque is 

C11TM ii 

2E.-~2 - - = E curl H, 

and the translational activity is 

— II curl E. 

Their sum is 

may, perhaps, be appropriately termed the Duplex method, since its characteristics are the exhibition of 
the electric, magnetic, and electromagnetic relations in a duplex form, symmetrical with respect to the 
electric and magnetic sides. But it is not merely a method of exhibiting the relations in a manner 
suitable to the subject, bringing to light useful relations which were formerly hidden from view by the 
intervention of the vector-potential and its parasites, but constitutes a method of working as well. 
There are considerable difficulties in the way of the practical employment of Maxwell's equations of 
propagation, even as they stand in his treatise. These difficulties are greatly magnified when we pro- 
ceed to more general cases, involving heterogeneity and eolotropy and motion of the medium supporting 
the fluxes. The duplex method supplies what is wanted. Potentials do not appear, at least initially. 
They are regarded strictly as auxiliary functions which do not represent any physical state of the 
medium. In special problems they may be of great service for calculating purposes; but in general 
investigations their avoidance simplifies matters greatly. The state of the field is settled by E and H, 
and these are the primary objects of attention in the duplex system. 

As the papers to which I have referred are not readily accessible, I may take this opportunity of 
mentioning that a Reprint of my { Electrical Papers ' is in the press (Macmillan and Co.), and that the 
first volume is nearly ready. 
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- aw VEH, 
making VEH the flux of energy.* 

All attempts to construct an elastic solid analogy with a distortional stress fail to 
give satisfactory results, because the energy is wrongly localised, and the flux of 
energy incorrect. Bearing this in mind, the above analogy is at first sight very 
enticing. But when we come to remember that the d/dt in (98) and (99) should be 
3/9 £, and find extraordinary difficulty in extending the analogy to include the conduc- 
tion current, and also remember that the electromagnetic stress has to be accounted 
for (in other words, the known mechanical forces), the perfection of the analogy, as 
far as it goes, becomes disheartening. It would further seem, from the explicit 
assumption that q = in obtaining W above, that no analogy of this kind can be 
sufficiently comprehensive to form the basis of a physical theory. We must go 
altogether beyond the elastic solid with the additional property of rotational elasticity. 
I should mention, to avoid misconception, that Sir W. Thomson does not push the 
analogy even so far as is done above, or give to /x and c the same interpretation. The 
particular meaning here given to [x is that assumed by Professor Lodge in his " Modern 
Views of Electricity," on the ordinary elastic solid theory, however. I have found it 
very convenient from its making the curl of the electric force be a Newtonian force 
(per unit volume). When impressed electric force e produces disturbances, their real 
source is, as I have shown, not the seat of e , but of curl e . So we may with facility 
translate problems in electromagnetic waves into elastic solid problems by taking the 
electromagnetic source to represent the mechanical source of motion, impressed New- 
tonian force. 

Examination of the Flux of Energy in a moving Medium, and Establishment of the 

Measure of " True " Current. 

§ !?• Now pass to the more general case of a moving medium with the equations 

curl H x = curl (H - h - h) = J = C + t + qp, (101) 

- curl E x = - curl (E - e - e) = G = K + B + ck (102) 

where E x is, for brevity, what the force E of the flux becomes after deducting the 
intrinsic and motional forces ; and similarly for Hp 
From these, in the same way as before, we deduce 

(e + e) J+ (ho + h) = EJ + HG + divVE^; ...... (103) 

and it would seem at first sight to be the same case again, but with impressed forces 

* This form of application of the rotating ether I ga,ve in ' The Electrician,' January 23, 1891, p. 360. 



446 ME. 0. HEAVISIDE 0E THE FORCES, STRESSES, AFD 

(e + e ) and (h + h ) instead of e and h 0? whilst the Poynting flux requires us to 
reckon only E x and H x as the effective electric and magnetic forces concerned in it.* 

But we must develop (Q + U + T) plainly first We have, by (86), (87), used in 
(103), 

e J + li G- = E (C -f- D + q P ) 4- H (K + B + q<x) - (eJ + hG) + div VE^. . . (104) 



Now here we have 



tr = y-|ED = J ED + i 

CvV 



~\ 



= ED + J (BE - ED) i 
= ED - * EiE 
= ED - TL. 



/^ • « 9 i o fl » 4 * V ^* ^ _/ 



_> 



* It will be observed that the constant 4tt, which usually appears in the electrical equations, is absent 
from the above investigations. This demands a few words of explanation. The units employed in the 
text are rational units, founded upon the principle of continuity in space of vector functions, and the 
corresponding appropriate measure of discontinuity, viz., by the amount of divergence. In popular 
language, the unit pole sends out one line of force, in the rational system, instead of 4w lines, as in the 
irrational system. The effect of the rationalisation is to introduce 4tt into the formula of central forces 
and potentials, and to abolish the swarm of 4tt's that appear in the practical formula of the practice of 
theory on Faraday-Maxwell lines, which receives its fullest and most appropriate expression in the 
rational method. The rational system was explained by me in i The Electrician/ in 1882, and applied to 
the general theory of potentials and connected functions in 1883. (Reprint, vol. 1, p. 199, and later, 
especially p. 262.) I then returned to irrational formulae because I did not think, then, that a reform of 
the units was practicable, partly on account of the labours of the JB. A. Committee on Electrical Units, 
and partly on account of the ignorance of, and indifference to, theoretical matters which prevailed at 
that time. But the circumstances have greatly changed, and I do think a change is now practicable. 
There has been great advance in the knowledge of the meaning of Maxwell's theory, and a diffusion of 
this knowledge, not merely amongst scientific men, but amongst a large body of practicians called into 
existence by the extension of the practical applications of electricity. Electricity is becoming, not only 
a master science, but also a very practical one. It is fitting, therefore, that learned traditions should not 
be allowed to control matters too greatly, and that the units should be rationalised. To make a begin- 
ning, I am employing rational units throughout in my work on. " Electromagnetic Theory," commenced in 
i The Electrician,' in January, 1891, and continued as fast as circumstances will permit ; to be republished 
in book form. In Section XVII. (October 16, 1891, p. 655), will be found stated more fully the nature 
of the change proposed, and the reasons for it. I p oint out, in conclusion, that as regards theoretical 
treatises and investigations, there is no difficulty in the way, since the connection of the rational 
and irrational units may be explained separately ; and I express the belief that when the merits of the 
rational system are fully recognised, there will arise a demand for the rationalisation of the practical 
units. We are, in the opinion of men qualified to judge, within a measurable distance of adopting 
the metric system in England. Surely the smaller reform I advocate should precede this. To put the 
matter plainly, the present system of units contains an absurdity running all through, it of the same 
nature as would exist in the metric system of common units were we to define the unit area to be 
the area of a circle of unit diameter. The absurdity is only different in being less obvious in the 
electrical case. It would not matter mnch if it were not that electricity is a practical science. 
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Comparison of the third with the second form of (105) defines the generalised mean- 
ing of 6 when c is not a mere scalar. Or thus, 

= KiE 1 8 + i<3 83 E a 8 + IW (106) 

representing the time-variation of U due to variation in the c's only. 
Similarly 

T == HB - £ H/iH == HB -■ T M > (107) 

with the equivalent meaning for //, generalised. 
Using these in (104) we have the result 

e J + h G = (Q + TJ + T) + q (E/> + H*) + (i E6E + | H/iH) - (eJ + MJ) + div VE^. . (108) 

Here we have, besides (Q + U + T), terms indicating the activity of a transla- 
tional force. Thus Ep is the force on electrification p> and Eqp its activity. Again, 



so that we have 



c 



u 



dc . 

dt 


+ qV.c ; 


do 

~ dt ' 


- qV.c, 

i 


dp 


1 



and, similarly, y (109) 



the generalised meaning of which is indicated by 



r\JT 

-^ + fEcE = - £E (qV . c) E = -- qVIT, ; (110) 



where, in terms of scalar products involving E and D, 



-qVTJ,= -|(E.qV.D-D.qV.E).. (m) 



This is also the activity of a translational force. Similarly. 



3t /x 

-^ + 1-HAH = - qVT„ . . . . (112) 



is the activity of a translational force. Then again 



- (eJ + h<j) = - JVqB - GVDq = q ( VJB + VDG) (113) 
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expresses a translational activity. Using them all in (108) it becomes 

e J + h G = (Q + U + T) + q(E/> + H<r - VTJ, - VT; + YJB + VDG) 

a 

+ div VE^ + g- (U, + T„,). .... (114) 

It is clear that we should make the factor of q be the complete translational force. 
But that has to be found; and it is equally clear that, although we appear to have 
exhausted all the terms at disposal, the factor of q in (114) is not the complete force, 
because there is no term by which the force on intrinsically magnetised or electrized 
matter could be exhibited. These involve e and h . But as we have 

q (Vj B + VDg ) = - (ej + hg ), .......... (115) 

a possible way of bringing them in is to add the left member and subtract the right 
member of (115) from the right member of (114) ; bringing the translational force to 
f, say, where 

f = E/> + H* - VU, - VT + V(J + j ) B + V (G + g ) D. • « • . • (116) 

But there is still the right number of (115) to be accounted for. We have 

- div (Velio + Ve h) = ej f tig + e^ + h g, ....... . (117) 

and, by using this in (114), through (115), (116), (117), we bring it to 

e J 4- h G = (Q + U+ T) + f q - (e j + h g) + cliv (VE^ - Veli - Ve li) + | (U c + T M ) ; (118) 



db 



or, transferring the e , h terms from the right to the left side. 



e J + h G = Q + U + T + fq + div (VE^ - Veh - Ve h) + K, (Pc + TJ . . (119) 

^ i vjlj 

Here we see that we have a correct form of activity equation, though it may not be 
the correct form. Another form, equally probable, is to be obtained by bringing in 
Veil ; thus 

diy Veil = h curl e — e curl h = — (ej + hg) = q (VjB + VDg) ? .... (120) 

which converts (119) to 

e J + h G n = Q + U + T + Fq + div (VE^ - Veh - Veh - Ve () h) + ^ (U c + TV) , (121) 

where F is the translational force 

F = E/3 + H<r- VU C - VT^ + Vcurl H.B + VcurlE.D, ... . (122) 
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which is perfectly symmetrical as regards E and H, and in the vector products 
utilises the fluxes and their complete forces, whereas former forms did this only 
partially. Observe, too, that we have only been able to bring the activity equation 
to a correct form (either (119) or (122)) by making e J be the activity of intrinsic 
force e , which requires that J should be the true electric current, according to the 
energy criterion, not J. 

§ 18. Now, to test (119) and (121), we must interpret the flux in (121), or say 



Y = VEiHi 



Veil - Veh - Ye li, 



« O ft • 



9 < O 



(123) 



which has replaced the Poynting flux VE^ when q = 0, along with the other 
changes. Since Y reduces to VE^ when q = 0, there must still be a Poynting 
flux when q is finite, though we do not know its precise form of expression. There 
is also the stress flux of energy and the flux of energy by convection, making a total 

X = W + q(U + T)~- 2% + q (U + T ), ....... (124) 



where W is the Poynting flux, and — 2Qg that of the stress, whilst q (U + T ) 
means convection of energy connected with the translational force. We should 
therefore have 



e J + h G = (Q + U + T) + (Q + lt + T ) + div X 
to express the continuity of energy. More explicitly 



e J + h G- = Q H- U + T + div [W + q (TJ + T)] 

+ Qo + U + T + div [~ 2% + q (IT + T )] 



(125) 



» » 



. (126) 



But here we may simplify by using the result (69) (with, however, f put = 0), making 
(126) become 

e J + h G = (Q + U + T) + Pq + Sa + div|> + q (U + T) - 2 Ctg], . . . (127) 

where S is the torque, and a the spin. 

Comparing this with (121), we see that we require 



W + q (U + T) ~~ 2% = VE^ - Veil - Ye li - Veh , 



9 * 



« • 



(128) 



with a similar equation when (119) Is used instead ; and we have now to separate the 
right member into two parts, one for the Poynting flux, the other for the stress flux, 
in such a way that the force due to the stress is the force F in (121), (122), or the 
force f in (119), (116); or similarly in other cases. It is unnecessary to give the 
failures; the only one that stands the test is (121), which satisfies it completely, 

MDCCOXOII. — A. 3 M 
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I argued that 



V {Mi " mm vnl \ XL "**""* Ma ) , » « « « » « » . » • ( JLiuy ) 



was the probable form, of the Poynting flux in the case of a moving medium, not 
VEjH^ because when a medium is endowed with a uniform translational motion, the 
transmission of disturbances through it takes place just as if it were at rest. With 
this expression (129) for W, we have, identically, 

VE^ - Veh - Yeoh - Velio == W - TeH - VEh. ...... (130) 

Therefore, by- (128) and (130), we get 

2ag = VeH + VEh + q(U + T), .......... (131) 

to represent the negative of the stress flux of energy, so that, finally, the fully 
significant equation of activity is 

e J + h G = Q + IJ + T + Fq + Sa + div [V (E - e ) (H - ii ) + q (U + T)] 

- cliv [VeH + VEh + q (U + T)]-. ........ (132) 

This is, of course, an identity, subject to the electromagnetic equations we started 
from, and is only one of the multitude of forms which may be given to it, many being 
far simpler. But the particular importance of this form arises from its being the 
only form apparently possible which shall exhibit the principle of continuity of energy 
without outstanding terms, and without loss of generality ; and this is only possible 
by taking J as the proper flux for e to work upon/" 

* In the original an erroneous estimate of the value of (dfdt>)(U c + IV) was used in some o£ the 
above equations. This is corrected. The following contains full details of the calculation. We require 
the value of (3/3*) TJ C , or of |E (dc/di)^ where do/dt is the linear operator whose components are the 
time- variations (for the same matter), of those of c. The calculation is very lengthy in terms of these 
six components. But veetorially it is not difficult. In (27), (28), we have 

D = cl = i.^E + j.CgE + k»e 3 E 1 q 

= (ix 1 4 j.c 2 + k.c 3 )E, j 
if the vectors C 13 c 2 , % are given by 

Cj = ic u + jc 13 + kc 13 , c 3 = ic gl + j% + k%» % = ic 31 + j%>, + ke 3B . 

¥e, therefore, have 

E S? -E = E( r -c +?.Co + ^^E + Bfi.^ + j.^+k.^)E. . . . (1325) 

dt \dt l dt H df V \ dt J dt dbj v ; 

The part played by the dots is to clearly separate the scalar products, 

Now snppose that the eolotropic property symbolised by o is intrinsically unchanged by the shift of 
the matter. The mere translation does not, therefore, affect it, nor does the distortion ; but the rotation 
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Derivation of the Electric and Magnetic Stresses and Forces from the Flux of 

Enerciv. 



§ 19. It will be observed that the convection of energy disappears by occurring 
twice oppositely signed ; but as it comes necessarily into the expression for the stress 
flux of energy, I have preserved the cancelling terms in (132). A comparison of the 
stress flux with the Poynting flux is interesting. Both are of the same form, viz., 
vector products of the electric and magnetic forces with convection terms ; but 
whereas in the latter the forces in the vector product are those of the field (i.e., only 
intrinsic forces deducted from E and H), in the former we have the motional forces 
e and h combined with the complete E and H of the fluxes. Thus the stress depends 



does. For if we turn round an eolotropic portion of matter, keeping* E unchanged, the value of U is 
altered by the rotation of the principal axes of c along with the matter, so that a torque is required. 

In equation (132a), then, to produce (132&), we keep E constant, and let the six vectors, i, j, k, c 1? c 3 , C 3 
rotate as a rigid body with the spin a = J curl q. Bat when a vector magnitude i is turned round in 
this way, its rate of time-change dijdt is Yai. Thus, for 9/8^, we may put Ya throughout. Therefore, 
by (1326), 



ic 



E ~~E = ECVaLCjL + Yaj.c 3 + Yak.c 3 )E + E (i.VaC]_ 4- j.Yac 3 -f k.Yac 3 )E. 

In this use the parallelepidedal transformation (12), and it becomes 

dc 
E ~- E = VEa(i.c 1 + j.c 3 + k.c 3 )E + E (i.c 1 -f j.c 3 + k,c 3 )YEa 

= (YEa)cE + Ec(YEa) = (D + D')YEa, . 



(132c) 



(lS2d) 



by (132a), if D' is conjugate to D ; that is, D' = c'E = Eo. 
have 



So, when c = c, as in the electrical case, we 



and similarly 



Pitt P\ 

V = i E 4 E = B YEa = aYDE, 

dt - at 

l L± _ i H /' H _ BVHa _ a VBH. 



>■ (132e) 



dt 



dt 



Now the torque arising from the stress is (see (139)) 



so we have 



3 
dt 



S = YDE + YBH, 



(U tf -f T^) = Sa = torque x spin. 



(132/) 



The variation allowed to i, j, k may seem to conflict with their constancy (as reference vectors) in 
general, But they merely vary for a temporary purpose, being fixed in the matter instead of in space. 
But we may, perhaps better, discard i, j, k altogether, and use any independent vectors, 1, m, n instead, 



making 



D = (I.Cjl + m.c 3 + n.c 3 ) E, 



• * 



• • (132gr) 



wherein the c's are properly chosen to suit the new axes. The calculation then proceeds as before, half 

8 M 2 
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entirely on the fluxes, however they be produced, in this respect resembling the 
electric and magnetic energies. 

To exhibit the stress, we have (131), or 

%i + daft + Qs& = VeH + YEh + q (IT + T) . (133) 

In this use the expressions for e and h, giving 

2% = YHYBq + YEYBq + q(U h T) 

= B.Hq - q.HB + D.Eq - q.EB + q (U + T) 

= (B.Hq - qT) + (D.Eq - qU) ; . . ........... . (134) 

where observe the singularity that q (IT + T) has changed its sign. The first set 
belongs to the magnetic, the second to the electric stress, since we see that the 
complete stress is thus divisible. 

The divergence of 2Qg being the activity of the stress-variation per unit volume, 
its N component is the activity of the stress per unit surface, that is, 

(NB.Hq - Kq.T) + (ND.Eq - Nq.TT) = q (H.BN + E.DN - NtJ - FT) = P N q. . (135) 
The stress itself is therefore 



the value of 3lT c /3£ arising from the variation of 1, m, n, and the other half from the C 3 s, provided c is 
irrotational. 

Or we may choose the three principal axes of c in the body, when 1, m, n will coincide with, and 
therefore move with them. 

Lastly, we may proceed thus :— 

3c 3d 3f 

E ^ E = E ^ - D ^- = EVaD -- DVaE = 2aVDE. ..... (1327.) 

That is, replace 3/3^ by Ya when the operands are E and D. This is the correct result, but it is not 
easy to justify the process directly and plainly; although the clue is given by observing that what we do 
is to take a difference, from which the time- variation of E disappears. 

If it is D that is kept constant, the result is 2aYED, the negative of the above. 

It is also worth noticing that if we split up E into "E x + E 3 we shall have 



E 1 |E a = a[V(E 1 e)E 2 ---VE 1 ( c E 2 )3 ; 



do 



*> ....... (132i) 



These are only equal when e = c', or Ec = cE ; so that, in the expansion of the torque, 

YDE = VDxEx + YD 3 E 3 + YD 3 E X + YD^, 

the cross-torques are not YDoE x and YD^, which are unequal, but are each equal to half the sum of 

these vector-products. 
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P N = (E.DN - NU) + (H.BN-NT), ......... (136) 

divided into electric and magnetic portions. This is with restriction to symmetrical 
[jb aad c, and with persistence of their forms as a particle moves, but is otherwise 
unrestricted. 

Neither stress is of the symmetrical or irrotational type in case of eolotropy, and 
there appears to be no getting an irrotational stress save by arbitrary assumptions 
which destroy the validity of the stress as a correct deduction from the electro- 
magnetic equations. But, in case of isotropy, with consequent directional identity 
of E and B, and of H and B, we see, by taking N in turns parallel to, or 
perpendicular to E in the electric case, and to H in the magnetic case, that the 
electric stress consists of a tension U parallel to E combined with an equal lateral 
pressure, whilst the magnetic stress consists of a tension T parallel to H combined 
with an equal lateral pressure. There are, in fact, Maxwell's stresses in an isotropic 
medium homogeneous as regards /x and c. The difference from Maxwell arises when 
/x and c are variable (including abrupt changes from one value to another of /x and c), 
and when there is intrinsic magnetisation, Maxwell's stresses and forces being then 
different. 

The stress on the plane whose normal is VEH, is 

E.DYEH + H.BVEH - (U + T) VEH 



Y EH 



E.HVDE + H.EVHB - (U -^ VEH 

, • • • • • » • • \ L*J t J 



V EH 



reducing simply to a pressure (U + T) in lines parallel to VEH in case of isotropy. 
§ 20. To find the force F, we have 

EN = div a N = div (D.EN - NIT + B.HN - NT) 
= EN./> + DV.EN - fE.NV.D -~ p.NV.E + &c 
= EN.p + D (V.EN - NV.E) + I (D.NV.E - E.NV.D) + &c. 
= N [E/> + V curl E.D - VU C + &c], (138) 

where the unwritten terms are the similar magnetic terms. This being the N 
component of F, the force itself is given by (122), as is necessary. j 

It is V curl h . B that expresses the translational force on intrinsically magnetised 
matter, and this harmonises with the fact that the flux B due to any impressed 
force h depends solely upon curl h . 

Also, it is ■— VT^ that explains the forcive on elastically magnetised matter, e.g., 
Fabaday's motion of matter to or away from the places of greatest intensity of the 
field, independent of its direction. 
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If S be the torque, it is given by 

VSN = P N - a N = E.DN - D.EM + &e. 

VN (VED + VHB) ; 

therefore 

S = YDE + VEH ............. (139) 

But the matter is put more plainly by considering the convergence of the stress- 
flux of enemy and dividing; it into traoslational and other parts. Thus 

div20g = Pq + (E.DV.q-Udivq) + (H.BV.q-Tdivq), (140) 

where the terms following Pq express the sum of the distortional and rotational 
activities. 

Shorter Way of going from the Circuital Equations to the Flux of Energy, Stresses, 

and Forces. 

§ 21. I have given the investigation in §§ 17 to 19 in the form in which it occurred 
to me before I knew the precise nature of the results, being uncertain as regards the 
true measure of current, the proper form of the Poynting flux, and how it worked in 
harmony with the stress flux of energy. But knowing the results, a short demonstra- 
tion may be easily drawn up, though the former course is the most instructive. Thus, 

start now from 

curl (H — ho) = J , 



• <J * » « 



(141) 

curl (E — e ) =: Gr , 

on the understanding that J and G are the currents which make e J and h Gr 
the activities of e and h the intrinsic forces. Then 

e J + h G = E J + HG + div W 5 (142) 

where 

W = V(E~~e )(H-h ); ........... (143) 

and we now take this to be the proper form of the Poynting flux, Now develop 
EJ and HG thus :— 

EJ + HG = E (C + D + qp + curl h) + H (K + B + q«r - GUr [ e ^ by ( 93 ) ; 

= Qi + U J r U c + Eq/o + E curl VDq 

+ Q2 + T -{- T^ + Hq<r + H curl VBq ? by (88) and (91) ; 

= Qi + XJ + to + Eq/> + E (D div q + qV.B - q div D - DV.q) 

+ Q 3 -f t + T^ + Hqo- + H (B div q + qV.B - q div B - BV.q), by (26), 

= Qi + U + U c + 2U div q + E.qV.D - E.DV.q 
-f magnetic terms, 

= (Qi + U + div qU) + (U div q - E.DV.q) + (jj c - qV.U + E.qV.D) 

4- magnetic terms. (144) 
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Now here 

qV.U = pl.qV.D + |D.qV.E, 

so that the terms in the third pair of brackets in (144) represent 

U 6 , + qV.U c = a 5 c = iE^E, 

ct ot 

with the generalised meaning before explained. So finally 

EJ + HO = Q + U + T + div q(U + T) + ^ (U, + T^) 

+ (U div q - E.DV.q) + (T div q - H.BV.q), . . . (145) 
which brings (142) to 

e J + h G = Q + IT + T + div {W + q (IT + T)} 

a 

+ *- (U c + T „) + (U div. q - E.DV.q) + (T div q - H.BV.q), . . . (146) 

which has to be interpreted in accordance with the principle of continuity of energy. 
Use the form (127), first, however, eliminating Fq by means of 

div2ttg = Eq +• 2aVq, 

which brings (127) to 

e J + h G = a + U + T + div {w + q (U + T)} - 2ttVg + Sa ; (147) 

and now, by comparison of (147) with (146) we see that 

3il 



Sa -f 2aVg = (E.DV.q - U div q) 



dt 



3T 

+ (H.BV.q - T div q) - -^ 5 ( 148 ) 

from which, when /x and do not change intrinsically, we conclude that 



a N = B.HN - NT + D.EN - NU 
P N = H.BN - NT + E.DN - NIT 






. (149) 



as before. In this method we lose sight altogether of the translational force which 
formed so prominent an object in the former method as a guide. 

Some Remarks on Hertz's investigation relating to the Stresses. 

§ 22. Variations of c and /x in the same portion of matter may occur in different 
ways, and altogether independently of the strain variations. Equation (146) shows 
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how their influence affects the energy transformations ; but if we consider only such 
changes as depend on the strain, i.e., the small changes of value which //, and c 
undergo as the strain changes, we may express them by thirty-six new coefficients 
each (there being six distortion elements, and six elements in /x, and six in c), and so 
reduce the expressions for d~U c /dt and dTJdi in (148) to the form suitable for 
exhibiting the corresponding change in Q N and in the stress function P N . As is 
usual in such cases of secondary corrections, the magnitude of the resulting formula 
is out of all proportion to the importance of the correction terms in relation to the 
primary formula to which they are added. 

Professor H. Hertz"' has considered this question, and also refers to von Helmoltz's 
previous investigation relating to a fluid. The c and p can then only depend on the 
density^ or on the compression, so that a single coefficient takes the place of the 
thirty-six. But I cannot quite follow Hertz's stress investigation. First, I would 
remark that in developing the expression for the distortional (plus rotational) activity, 
he assumes that all the coefficients of the spin vanish identically ; this is done in 
order to make the stress be of the irrotational type. But it may easily be seen that 
the assumption is inadmissible by examining its consequence, for which we need only 
take the case of c and /a intrinsically constant- By (139) we see that it makes S = 0, 
and therefore (since the electric and magnetic stress are separable), VHB = 0, and 
VED = ; that is, it produces directional identity of the force E and the flux D, 
and of the force H and the flux B. This means isotropy, and, therefore, breaks 
down the investigation so far as the eolotropic application, with six fi and six c 
coefficients, goes. Abolish the assumption made, and the stress will become that used 
by me above. 

Another point deserving of close attention in Hertz's investigation, relates to the 
principle to be followed in deducing the stress from the electromagnetic equations. 
Translating into my notation it would appear to amount to this, the a priori assump- 
tion that the quantity 

— ~^r l ^ n s *■ * * • ' ' * * * • • ♦ ( LOU ) 

v dt\ / . . J 

where v indicates the volume of a moving unit element undergoing distortion, may 
be taken to represent the distortional (plus rotational) activity of the magnetic 
stress. Similarly as regards the electric stress. 
Expanding (150) we obtain 

v-l-^rrETT-hT cllV Q — V^ ......... (151) 

dt v dt dt H dt y J 

Now the second circuital law (90) may be written 



- curl (E - e ) = K + y + (B div q - BV.q)- . • • . . . . • (152) 
# ' Wiedemann's Annalen,' v. 41, p. 869, 
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Here ignore e , K, and ignore the curl of the electric force, and we obtain, by using 
(152) in (151), 

H.BV.q - HB div q + T div q - —^ = H.BV.q - T div q ~ ^, • • • (153) 

which represents the distortional activity (my form, not equating to zero the 
coefficients of curl q in its development). We can, therefore, derive the magnetic 
stress in the manner indicated, that is, from (150), with the special meaning of dB/dt 
later stated, and the ignorations or nullifications, 

In a similar manner, from the first circuital law (89), which may be written 

curl (H - h ) = C + y + (D div q - DV.q), ....... (154) 

we can, by ignoring the conduction current and the curl of the magnetic force, obtain 

l|(,U )== E.DV. q -Udi vq -^, (155) 

which represents the distortional activity of the electric stress. 

The difficulty here seems to me to make it evident & priori that (150), with the 
special reckoning of dB/dt should represent the distortional activity (plus rotational 
understood) ; this interesting property should, perhaps, rather be derived from the 
magnetic stress when obtained by a safe method. The same remark applies to the 
electric stress. Also, in (150) to (155) we overlook the Poyntjng flux. I am not 
sure how far this is intentional on Professor Hertz's part, but its neglect does not 
seem to give a sufficiently comprehensive view of the subject. 

The complete expansion of the magnetic distortional activity is, in fact, 

3t 

H.BV.q - T div q - ^ = Q 3 + T + div qT - HG ; (156) 

and similarly, that of the electric stress is 

E.DV.q - U div q — ^ = Qi + U + div qU - EJ (157) 

It is the last term of (156) and the last term of (157), together, which bring in the 
Poynting flux. Thus, adding these equations, 

2ttV 2 - jk(U e + T„) = Q + U + t + div q (IT + T) - (EJ + HG ), . . . (158) 

where 

(E J + HGr„) = (e„J + h Gr ) — div W ; (1 59) 

MDCCCXCII. — A. 3 N 
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and so we come round to the equation of activity again, in the form (146), by using 
(159) in (158). 

Modified Form of Stress-vector, and Application to the Surface separating tivo Regions. 

§ 23* The electromagnetic stress, P N of (149) and (136) may be put into another 
interesting form. We may write it 

P N = i (E.ND + V.VHE.B) + | (H.NB 4- Y.YHH.E) (160) 

Now, NX) is the surface equivalent of div D and ¥B of div B ; whilst VNE and VNH 
are the surface equivalents of curl E and curl IL We may, therefore, write 



^ = i (E P ' + YDGT) + | (HV + VJ'B), ........ (161) 



and this is the force, reckoned as a pull, on unit area of the surface whose normal is 
N- Here the accented letters are the surface equivalents of the same quantities 
unaccented, which have reference to unit volume. 

Comparing with (122) we see that the type is preserved, except as regards the 
terms in F due to variation of c and /x, in space. That is, the stress is represented in 
(101) as the translational force, due to E and H, on the fictitious electrification, 
magnetification, electric current, and magnetic current produced by imagining E and 
H to terminate at the surface across which P N is the stress. 

The coefficient | which occurs in (161) is understandable by supposing the fictitious 
quantities (" matter " and " current ") to be distributed uniformly within a very thin 
layer, so that the forces E and H which act upon them do not then terminate quite 
abruptly, but fall off gradually through the layer from their full values on one side to 
zero on the other. The mean values of E and H through the layer, that is, ^E and 
-|H are thus the effective electric and magnetic forces on the layer as a whole, per 
unit volume density of matter or current ; or ^ and -pi per unit surface density 
when the layer is indefinitely reduced in thickness. 

Considering the electric field only, the quantities concerned are electrification and 
magnetic current. In the magnetic field only they are magnetification and electric 
current. Imagine the medium divided into two regions A and B, of which A is 
internal, B external, and let N be the unit normal from the surface into the external 
region. The mechanical action between the two regions is fully represented by the 
stress P N over their interface, and the forcive of B upon A is fully represented by 
the E and H in B acting upon the fictitious matter and current produced on the 
boundary of B, on the assumption that E and H terminate there. If the normal 
and P N be drawn the other way, thus negativing them both, as well as the fictitious 
matter and current on the interface, then it is the forcive of A on B that is repre- 
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sented by the action of E and H in A on the new interfacial matter and current. 
That is, the E and H in the region A may be done away with altogether, because 
their abolition will immediately introduce the fictitious matter and current equivalent, 
so far as B is concerned, to the influence of the region A. Similarly E and H in B 
may be abolished without altering them in A. And, generally, any portion of the 
medium may be taken by itself and regarded as being subjected to an equilibrating 
system of forces, when treated as a rigid body. 




§ 24. When c and jjl do not vary in space, we do away with the forces -— -|E 2 Vc 
and -— ^H 2 V/i, and make the form of the surface and volume translational forces 
agree. We may then regard every element of p or of cr as a source sending out 
from itself displacement and induction isotropically, and every element of J or G- as 
causing induction or displacement according to Ampere's rule for electric current and 
its analogue for magnetic current. Thus 

E = 2 p/ ° ± Y J X ~ , • • • (162) 

H = 2 ^ + VJr 1; 

where r l is a unit vector drawn from the infinitesimal unit volume in the summation 
to the point at distance r where E or H is reckoned. Or, introducing potentials, 

(164) 



E = ■ 


-V2? /C -curl 2, G 
4?7T 4nrr 


11 = 


- V2 ^ + curl 2 / 

4s7TV 4&TTT 



(165) 



These apply to the whole medium, or to any portion of the same, with, in the 
latter case, the surface matter and current included, there being no E or I outside 
the region, whilst within it E and H are the same as due to the matter and current 
in the whole region ("matter," p and a ; "current," J and G). But there is no 
known general method of finding the potentials when c and p vary. 

We may also divide E and H into two parts each, say E x and H x due to matter 
and current in the region A, and E 2 , H 2 due to matter and current in the region B 

3 N 2 
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surrounding it, determinable in the isotropic homogeneous case by the above formulae. 
Then we may ignore E 2 and H x in estimating the forcive on the matter and current 
in the region A ; thus, 

(Hg^ + VJiBsj) +2(E^ 1 .+ VD 2 G 1 ), . ........ (166) 



where a l = div B l = div B, and J x = curl Hj = curl H in region A, is the resultant 
force on the region A, and 

2(1^3 + VJgB JL ) + 2(E 1 ^ + VDiGg), . . . (167) 

is the resultant force on the region B ; the resultant force on A due to its own E and 
H being zero, and similarly for B. These resultant forces are equal and opposite, and 
so are the equivalent surface-integrals 

2(H 3 «r 1 ' + VJ/B 3 )+2(E ?J > 1 '+YD 3 G 1 '), ........ (168) 

and 

2 (H lff8 ; + VJ 8 'B X ) +2(E l/ > 3 / + VD 1 G 3 '), ........ (169) 

taken over the interface. The quantity summed is that part of the stress- vector, P N , 
which depends upon products of the H of one region and the B of the other, &c. 
Thus, for the magnetic stress only, 

H.BN - N.iHB = (H^H - H.iH^) + (H^N - N.^HA) 

+ (Hsj.B 2 N - N.|H 3 B 3 ) + (H^N - N^H^), . . . (170) 

and it is the terms in the second and fourth brackets (which, be it observed, are not 
equal) which together make up the magnetic part of (168) and (169) or their nega- 
tives, according to the direction taken for the normal ; that is, since H 1 B 2 = H^B l5 

2P N = 2(H 1 .B 3 H + E^N - MA) = 2(h.BN - N.pffl) 

= 2(H 1 <t 2 / + VJjj'Bi) = 2(h 3 v + VJ/B 2 ) = 2(bV + VJ'B) 
==2F = 2(H 1 (r 2 + VJ5 J B 1 )==2(H 30 - 1 + VJ l B 2 ) = 2(H < r + VJB), . . . (171) 

where the first six expressions are interfacial summations, and the four last summa- 
tions throughout one or the other region, the last summation applying to either 
region. No special reckoning of the sign to be prefixed has been made. The 
notation is such that H = H x + H 2 , or = <r l + cr 2 , &c, &c. 

The comparison of the two aspects of electromagnetic theory is exceedingly curious ; 
namely, the precise mathematical equivalence of " explanation " by means of instan- 
taneous action at a distance between the different elements of matter and current, 
each according to its kind, and by propagation through a medium in time at a finite 
velocity. But the day has gone by for any serious consideration of the former view 
other than as a mathematical curiosity, 
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Quaternionic Form of Stress- Vector. 

t , ^u. We may also notice the Quaternion form for the stress function, which is so 
vital a part of the mathematics of forces varying as the inverse square of the distance, 
and of potential theory. Tsotropy being understood, the electric stress may be 
written 

P N = ic[EN" 1 E], (172) 

where the quantity in the square brackets is to be understood quaternionically. It 
is, however, a pure vector. Or, 



N 



E 



c 

2 



(173) 



E 



that is, not counting the factor \ c, the quaternion ~ is the same as the quaternion 

; or the same operation which turns "N to E also turns B to P N . Thus, N, E, 

and P N are in the same plane, and the angle between N and E equals that between 
E and P N ; and E and P N are on the same side of N when E makes an acute angle 
with N. Also, the tensor of P N is IT, so that its normal and tangential components 

A 

are U cos 20 and U sin* 20, if = NE. 
Otherwise 

N =— ^<;[ENE], ............ (174) 

since the quaternionic reciprocal of a vector has the reverse direction. The corre- 
sponding volume translational force is 

F = - cV [EVE], (175) 

which is also to be understood quaternionically, and expanded, and separated into 
parts to become physically significant. I only use the square brackets in this 
paragraph to emphasise the difference in notation, It rarely occurs that any 
advantage is gained by the use of the quaternion, in saying which, I merely repeat 
what Professor Willakd Gibbs has been lately telling us ; and I further believe the 
disadvantages usually far outweigh the advantages. Nevertheless, apart from practical 
application, and looking at it from the purely quaternionic point of view, I ought to 
also add that the invention of quaternions must be regarded as a most remarkable 
feat of human ingenuity. Vector analysis, without quaternions, could have been 
found by any mathematician by carefully examining the mechanics of the Cartesian 
mathematics ; but to find out quaternions required a genius. 
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Remarks on the Translational Force in Free Ether. 

§ 26„ The little vector Veil, which has an important influence in the activity 
equation, where e and h are the motional forces 

e = VqB, h = VDq, 

has an interesting form, viz., by expansion, 

Veil = q.qVDB-= \ • qVEH, . . . . (176) 

if v be the speed of propagation of disturbances. We also have, in connection there- 
with, the equivalence 

always. 

The translational force in a non-conducting dielectric, free from electrification and 

intrinsic force, is 

F = VJB + YDG + VjB + VDg, 
or, approximately, 



= VDB + VDB = ~ VDB == 4 4 VEH = Z1 9 (178) 

at v" at v^ v y 

The vector VDB, or the flux of energy divided by the square of the speed of 
propagation, is, therefore, the momentum (translational, not magnetic, which is quite 
a different thing), provided the force P is the complete force from all causes acting, 
and we neglect the small terms VjB and VDg 8 . 

But have we any right to safely write 



F = m^~, . . . . . . . . . . . . (179) 

where m is the density of the ether ? To do so is to assume that F is the only force 
acting, and, therefore, equivalent to the time-variation of the momentum of a moving 
particle. # 

Now, if we say that there is a certain forcive upon a conductor supporting electric 
current ; or, equivalently, that there is a certain distribution of stress, the magnetic 
stress, acting upon the same, we do not at all mean that the accelerations of momentum 
of the different parts are represented by the translational force, the "electromagnetic 
force." It is, on the other hand, a dynamical problem in which the electromagnetic 
force plays the part of an impressed force, and similarly as regards the magnetic 

# Professor J. J. Thomson has endeavoured to make practical use of the idea, 'Phil. Mag./ March, 
1891. See also my article, ' The Electrician,' January 15, 1886, 
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stress ; the actual forces and stresses being only determinable from a knowledge of 
the mechanical conditions of the conductor, as its density, elastic constants, and the 
way it is constrained. Now, if there is any dynamical meaning at all in the electro- 
magnetic equations, we must treat the ether in precisely the same way. But we do 
not know, and have not formularised, the equations of motion of the ether, but only 
the way it propagates disturbance through itself, with due allowance made for the 
effect thereon of given motions, and with formularisation of the reaction between the 
electromagnetic effects and the motion. Thus the theory of the stresses and forces in 
the ether and its motions is an unsolved problem, only a portion of it being known so 
far, i.e., assuming that the Maxwellian equations do express the known part. 

When we assume the ether to be motionless, there is a partial similarity to the 
theory of the propagation of vibrations of infinitely small range in elastic bodies, when 
the effect thereon of the actual translation of the matter is neglected. 

But in ordinary electromagnetic phenomena, it does not seem that the ignoration 
of q can make any sensible difference, because the speed of propagation of disturbances 
through the ether is so enormous, that if the ether were stirred about round a magnet, 
for example, there would be an almost instantaneous adjustment of the magnetic 
induction to what it would be were the ether at rest. 

Static Consideration of the Stresses. — Indeterminateness. 

§ 27. In the following the stresses are considered from the static point of view, 
principally to examine the results produced by changing the form of the stress func- 
tion. Either the electric or the magnetic stress alone may be taken in hand. Start 
then, from a knowledge that the force on a magnetic pole of strength m is Em, where 
E. is the polar force of any distribution of intrinsic magnetisation in a medium, the 
whole of which has unit inductivity, so that 

div R = m = conv h (180) 

measures the density of the fictitious "magnetic" matter; h being the intrinsic 
force, or, since here /x = 1, the intensity of magnetisation. The induction is 
B = h + B. This rudimentary theory locates the force on a magnet at its poles, 

superficial or internal, by 

E = E div E (181) 

The N component of F is 

FN = EN.divE = div{E.EN-N.|E 3 }, (182) 

because curl K = 0. Therefore 

P N = E.EN -~ N.^E 3 (183) 

is the appropriate stress, of irrotational type. Now, however uncertain we may be 
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about the stress in the interior of a magnet, there can be no question as to the 
possible validity of this stress in the air outside our magnet, for we know that the 
force E is then a polar force, and that is all that is wanted, m and h being merely 
auxiliaries, derived from E. 

Now consider a region A, containing magnets of this kind, enclosed in B, the rest 
of space, also containing magnets. The mutual force between the two regions is 
expressed by 2P N over the interface, which we may exchange for 2 Em through 
either region A or B, still on the assumption that E remains polar. 

But if we remove this restriction upon the nature of E, and allow it to be arbitrary, 
say in region B or in any portion thereof, we find 

HE = div P N = EH div R + NV curl R . E; 
or 

F == Em + YJR, 

if J = curl R. This gives us, from a knowledge of the external magnetic field of polar 
magnets only, the mechanical force exerted by a magnet on a region containing J, 
whatever that may be, provided it be measurable as above ; and without any experi- 
mental knowledge of electric currents, we could now predict their mechanical effects in 
every respect by the principle of the equality of action and reaction, not merely as 
regards the mutual influence of a magnet and a closed current, but as regards the 
mutual influence of the closed currents themselves; the magnetic force of a closed current, 
for instance, being the force on unit of m, is equivalently the force exerted by m on the 
closed current, which, by the above, we know. Also, we see that according to this 
magnetic notion of electric current, it is necessarily circuital. 

At the same time, it is to be remarked that our real knowledge must cease at the 
boundary of the region containing electric current, a metallic conductor for instance ; 
the surface over which P N is reckoned, on one side of which is the magnet, on the 
other side electric current, can only be pushed up as far as the conductor. The stress 
P N may therefore cease altogether on reaching the conductor, where it forms a distri- 
bution of surface force fully representing the action of the magnet on the conductor. 
Similarly, we need not continue the stress into the interior of the magnet. Then, so 
far as the resultant force on the magnet as a whole, in translating or rotating it, and, 
similarly, so far as the action on the conductor, is concerned, the simple stress P N of 
constant tensor -|E 3 , varying from a tension parallel to E to an equal pressure laterally, 
acting in the medium between the magnet and conductor, accounts, by its terminal 
pulls or pushes, for the mechanical forces on them. The lateral pressure is especially 
prominent in the case of conductors, whilst the tension goes more or less out of sight, 
as the immediate cause of motion. Thus, when parallel currents appear to attract one 
another, the conductors are really pushed together by the lateral pressure on each 
conductor being greater on the side remote from the other than on the near side : 
whilst if the currents are oppositely directed, the pressure on the near sides is greater 
than on the remote sides, and they appear to repel one another, 
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The effect of continuing the stress into the interior of a conductor of unit induc- 
tivity, according to the same law, instead of stopping it on its boundary, is to 
distribute the translational force bodily, according to the formula 2VJK, instead of 
superficially, according to 2P N . In either case, of course, the conductor must be 
strained by the magnetic stress, with the consequent production of a mechanical stress. 
But the strain (and associated stress) will be different in the two cases, the applied 
forces being differently localised. The effect of the stress on a straight portion of a 
wire supporting current, due to its own field only, is to compress it laterally, and to 
lengthen it. Besides this, there will be resultant force on it arising from the different 
pressures on its opposite sides due to the proximity of the return conductor or rest of 
the circuit, tending to move it so as to increase the induction through the circuit per 
unit current, that is, the inductance of the circuit, 

§ 28. If now, we bring an elastically magnetisable body into a magnetic field, it 
modifies the field by its presence, causing more or less induction to go through it 
than passed previously in the air it replaces, according as its inductivity exceeds or is 
less than that of the air. The forcive on it, considered as a rigid body, is completely 
accounted for by the simple stress P N in the air outside it, reckoned according to the 
changed field, and supposed to terminate on the surface of the disturbing body. This 
is true whether the body be isotropic or heterotopic in its inductivity ; nor need the 
induction be a linear function of the magnetic force. It is also true when the body is 
intrinsically magnetised ; or is the seat of electric current. In short, since the 
external stress depends upon the magnetic force outside the body, when we take the 
external field as we may find it, that is, as modified by any known or unknown causes 
within the body, the corresponding stress, terminated upon its boundary, fully 
represents the forcive on the body, as a whole, due to magnetic causes. This follows 
from the equality of action and reaction ; the force on the body due to a unit pole is 
the opposite of that of the body on the pole. 

If we wish to continue the stress into the interior of the body, surrounded on all 
sides by the unmagnetised medium of unit inductivity, as we must do if we wish to 
arrive ultimately at the mutual actions of its different parts, and how they are modified 
by variations of inductivity, by intrinsic magnetisation,, and by electric current in the 
body, we may, so far as the resultant force and torque on it are concerned, do it in 
any way we please, provided we do not interfere with the stress outside. For the 
internal stress, of any type, will have no resultant force or torque on the body, and 
there is merely left the real external stress. 

Practically, however, we should be guided by the known relations of magnetic 
force, induction, magnetisation, and current, and not go to work in a fanciful manner; 
furthermore, we should always choose the stress in such a way that if, in its expres- 
sion, we take the inductivity to be unity, and the intrinsic magnetisation zero, it 
must reduce to the simple Maxwellian stress in air (assumed to represent ether here). 

mbcccxcii. — A. 3 o 
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But as we do nob know definitely the forcive arising from the magnetic stress in the 
interior of a magnet, there are several formulae that suggest themselves as possible. 

Special Kinds of Stress Formulae statically suggested. 

§ 29. Thus, first we have the stress (183) ; let this be quite general, then 

fP N = RitN-N,|R2, (184) 

[F = RdivR + VJR (185) 

Here R is the magnetic force of the field, not of the flux B. If /x = 1, div R is 
the density of magnetic matter, the convergence of the intrinsic magnetisation, but 
not otherwise. In general, it is the density of the matter of the magnetic potential, 
calculated on the assumption /x = 1. The force on a magnet is located in this system 
at its poles, whether the magnetisation be intrinsic or induced. The second term in 
(185) represents the force on matter bearing electric current (J = curl R), but has to 
be supplemented by the first term, unless div R = at the place. 

§ 30. Next, let the stress be /x times as great for the same magnetic force, but be 
still of the same simple type, /x being the inductivity, which is unity outside the body, 
but having any positive value, which may be variable, within it. Then we shall have 

f P N = R.N//R - N.iR/iR, (186) 

(2) < 

[F = Rw 4- VJ/iR - i-R 2 V//, (187) 

where m = conv /xh = div /xR is the density of magnetic matter, /xh being the 
intensity of intrinsic magnetisation. 

The electromagnetic force is made /x times as great for the same magnetic force ; 
the force on an intrinsic magnet is at its poles ; and there is, in addition, a force 
wherever /x varies, including the intrinsic magnet, and not forgetting that a sudden 
change in /x, as at the boundary of a magnet, has to count. This force, the third term 
in (187), explains the force on inductively magnetised matter. It is in the direction 
of most rapid decrease of /x. 

§ 31. Thirdly, let the stress be of the same simple type, but taking H instead of 
R, H being the force of the flux B = /xH = /x (R + h ), where h is as before. 

We now have 

fP N = H.NB - N.iHB, (188) 

LP = VJB + Vj () B - |H 2 V,<, (180) 

where j = curl h is the distribution of fictitious electric current which produces the 
same induction as the intrinsic magnetisation /xh , and J is, as before, the real current. 
It is now gwasi-electromagnetic force that acts on an intrinsic magnet, with, 
however, the force due to V/a, since a magnet has usually large /x compared with air. 
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The above three stresses are all of the simple type (equal tension and perpendicular 
pressure), and are irrotational, unless /x be the eolotropic operator. No change is, in 
the latter case, needed in (186), (188), whilst in the force formulae (187), (189), the 
only change needed is to give the generalised meaning to V/x. Thus, in (1 89), instead 
of H 2 V/x, use 2V1V, 

or 

V /x (H/iH), 

or 

i (»* , ) t '( , * , ) ti ( H W 

or 

(v„-v„)hb, 

or 

i (HViB - BV X H) + j (HVoB - BV.H) + k (HV 3 B - BV S H), 

showing the i, j, k components. 

Similarly in the other cases occurring later. 

The following stresses are not of the simple type, though all consist of a tension 
parallel to R or H combined with an isotropic pressure. 

§ 32. Alter the stress so as to locate the force on an intrinsic magnet bodily upon 
its magnetised elements. Add R./xh N to the stress (186), and therefore /xh () .RN 
to its conjugate ; then the divergence of the latter must be added to the N-component 
of the force (187). Thus we get, if I = /xh , 

f P N = R.BN - N,| R//R, (190) 

W < 

[F = IV.R + VJ/iR - -jLR*V/« 091) 

But here the sum of the first two terms in F may be put in a different form. Thus, 

IV.R = I^R + LVoR + I S V 3 R 
= i.IVRi + j.IVlt- 2 + k.IVR s . 

Also 

IVlli = IVjR + I (VR-! - V X R) = IViR + iVJI. 
These bring (191) to 

F = (i.IViR + j.IVoR + k.IV.,R) + VJB - i R 2 V//, (192) 

where the first component (the bracketted part) is Maxwell's force on intrinsic 
magnetisation, and the second his electromagnetic force. The third, as before, is 
required where //, varies. 

§ 33. To the stress (190) add — N.-^RJ, without altering the conjugate stress 
making 

3 o 2 
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\ P ~ VJB - i {i (RV X B - BVjR) + j (EV 2 B - BV 3 E) + k (RV 3 B - BV S R)} . . (194) 
= VJB-(V D -V E )pB. 

This we need not discuss, as it is merely a transition to the next form. 
§ 34. To the stress (193) add h .NB ; we then get 

J iTj^ ^^ Ji.JMJj ~~~ JM.i^HJd, ............. (±yo) 

(6) <j F = YJB + fi.BV/^ -f j.BVfe 3 + k.BVA 3 } 

L - | {i. (EV^ - BV X R) + j. (RV 3 B - BV 3 R) + k. (RV 8 B - BV 3 R)} 5 . . . . (196) 
= YJB + BV.h - (V B - V R ) |RB, 

where h i} h 2 , h z are the components of h . 

Now if to this last stress (195) we add — N.|h B, we shall come back to the 
third stress, (188), of the simple type. 

Perhaps the most instructive order in which to take the six stresses is (1), (2), (4), 
(5), (6), and (3) ; merely adding on to the force, in passing from one stress to the 
next, the new part which the alteration in the stress necessitates, 
. To the above we should add Maxwell's general stress, which is 

(7) J F == YJB + {iJViR + j.IV 2 R + kJV 3 R} 

[ + {LMTC^R + j.MV 3 R + k.MV 3 RJ, ...... (198) 

= YJB + V R [R(I+M)] ; 

where M = (/x, — l) R = intensity of induced magnetisation. There is a good deal 
to be said against this stress ; some of which later. 



Remarks on Maxwell's General Stress. 

§ 35. All the above force formulae refer to the unit volume ; whenever, therefore, 
a discontinuity in the stress occurs at a surface, the corresponding expression per unit 
surface is needed; i.e., in making a special application, for it is wasted labour else. 
It might be thought that as Maxwell gives the force (198), and in his treatise 
usually gives surface expressions separately, so none is required in the case of this his 
force system (198). But this formula will give entirely erroneous results if carried 
out literally. It forms no exception to the rule that all the expressions require 
surface additions. 

Maxwell's general stress has the apparent advantage of simplicity. It merely 
requires an alteration in the tension parallel to R, from R 2 to RB, whilst the lateral 
pressure remains |R 3 , when we pass from unmagnetised to magnetised matter. The 
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force to which it gives rise is also apparently simple, being merely the sum of two 
forces, one the electromagnetic, VJB, the other a force on magnetised matter whose i 
component is (I + M)V 1 R, both per unit volume, the latter being accompanied (in 
case of eolotropy) by a torque. Now I is the intrinsic and M the induced magneti- 
sation, so the force is made irrespective of the proportion in which the magnetisation 
exists as intrinsic or induced. In fact, Maxwell's " magnetisation " is the sum of 
the two without reservation or distinction. But to unite them is against the whole 
behaviour of induced and intrinsic magnetisation in the electromagnetic scheme of 
Maxwell, as I interpret it. Intrinsic magnetisation (using Sir W. Thomson's term) 
should be regarded as impressed (I = /xh , where h is the equivalent impressed 
magnetic force); on the other hand, " induced" magnetisation depends on the force 
of the field {M = (/x — 1)R}. Intrinsic magnetisation keeps up a field of force. 
Induced magnetisation is kept up by the field. In the circuital law I and M 
therefore behave differently. There may be absolutely no difference whatever 
between the magnetisation of a molecule of iron in the two cases of being in a 
permanent or a temporary magnet. That, however, is not in question. We have no 
concern with molecules in a theory which ignores molecules, and whose element of 
volume must be large enough to contain so many molecules as to swamp the charac- 
teristics of individuals. It is the resultant magnetisation of the whole assembly that 
is in question, and there is a great difference between its nature according as it 
disappears on removal of an external cause, or is intrinsic. The complete amalgama- 
tion of the two in Maxwell's formula must certainly, I think, be regarded as a false 
step. 

We may also argue thus against the probability of the formula. If we have a 
system of electric current in an unmagnetisable (/jl = 1) medium, and then change /x 
everywhere in the same ratio, we do not change the magnetic force at all, the 
induction is made /x times as great, and the magnetic energy /x times as great, and is 
similarly distributed. The mechanical forces are, therefore, /x times as great, and are 
similarly distributed. That is, the translational force in the /x = 1 medium, or VJR, 
becomes VJ/xR in the second case in which the inductivity is /x, without other 
change. But there is no force brought in on magnetised matter per se. 

Similarly, if in the /x = 1 medium we have intrinsic magnetisation I, and then 
alter /x in any ratio everywhere alike, keeping I unchanged, it is now the induction 
that remains unaltered, the magnetic force becoming /x""" 1 times, and the energy /x" 1 
times the former values, without alteration in distribution (referring to permanent 
states, of coui*se). Again, therefore, we see that there is no translational force 
brought in on magnetised matter merely because it is magnetised. 

Whatever formula, therefore, we should select for the stress function, it would 
certainly not be Maxwell's, for cumulative reasons. When, some six years ago, I 
had occasion to examine the subject of the stresses, I was unable to arrive at any 
very definite results, except outside of magnets or conductors. It was a perfectly 
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indeterminate problem to find the magnetic stress inside a body from the existence of 
a known, or highly probable, stress outside it. All one could do was to examine the 
consequences of assuming certain stresses, and to reject those which did not work 
well After going into considerable detail, the only two which seemed possible were 
the second and third above (those of equations (186) and (188) above). As regards 
the seventh (Maxwell's stress equation (198) above), the apparent simplicity produced 
by the union of intrinsic and induced magnetisation, turned out, when examined into 
its consequences, to lead to great complication and unnaturalness. This will be 
illustrated in the following example, a simple case in which we can readily and fully 
calculate all details by different methods, so as to be quite sure of the results we 
ought to obtain. 

A ivorked-out Example to Exhibit the Forcives contained in Different Stresses. 

§ 36. Given a fluid medium of inductivity /x l5 in which is an intrinsic magnet of 
the same inductivity. Calculate the attraction between the magnet and a large solid 
mass of different inductivity /x 3 . Here it is only needful to calculate the force on a 
single pole, so let the magnet be infinitely thin and long, with one pole of strength m 
at distance a from the medium /x 3 , which may have an infinitely extended plane 
boundary. By placing a fictitious pole of suitable strength at the optical image in 
the second medium of the real pole in the first, we may readily obtain the solution. 




Let PQ be the interface and the real pole be at A and its image at B. We have 
first to calculate the distribution of R, magnetic force, in both media due to the 
pole m, as disturbed by the change of inductivity. We have div /^Rj = m in the 
first medium, and div /x 2 R 2 — ' m ^e second, therefore R has divergence only on 
the interface. Let or be the surface density of the fictitious interfacial matter to 
correspond ; its force goes symmetrically both ways ; the continuity of the normal 
induction therefore gives, at distance -r from A, the condition 
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/'if-; — — ~< "" i") = ^\~a 1 + ^ '» (199) 

because m/47rjLt 1 r 2 is the tensor of the magnetic force due to m in the fi x medium when 
of infinite extent. Therefore 



(T 



__ /e l — /* 2 7WY.I 



(200) 



The magnetic potential XI, such that R = — Vfl is the polar force in either region, is 
therefore the potential of / m/fi i at A and of a over the interface. 
But if we put matter n at the image B, of amount 

n=/V^O! 1 , (201) 

/*1 + /*2 / l l 

the normal component of R 2 on the /x 2 side due to n and the pole rn will be 



ma na ma ^ 

4^1^ 47T? 3 49T/4J?" 3 



(202) 



the same value as before ; the force R x on the /^ side is, therefore, the same as that 
due to matter m//^ at A and matter n at B ; whilst on the /x, 2 side the force R 2 is 

2m 
that due to matter m/u, at A and matter n also at A, that is, to matter at A. 

Thus in the /x 2 medium the force R 2 is radial from A as if there were no change of 
indue tivity, though altered in intensity. 

The repulsion between the pole m and the solid mass is not the repulsion between 
the matters ra//^ and n of the potential, but is 

= m x magnetic force at A due to matter n at B, 
= n X magnetic force at B due to matter m//^ at A, 

~ 4>r (2a)* ^T^o* 47^(2^' ' • ^ 6) 

becoming an attraction when /x 2 > /jl 19 making n negative. When /x 2 = 0, the 
repulsion is 

4^" (2aY ' 



m 2 



when /a 2 = oo, it is turned into an attraction of equal amount. 

Similarly, if we consider the attraction to be the resultant force between m and the 
interfacial matter a, we shall get the same result by 



(Tina 



2^, (204) 

47rr* 



the quantity summed (over the interface) being <r X normal component of magnetic 
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force due to matter m in a medium of unit inductivity, or the normal component of 
induction due to m in its own medium. For this is 



i 



ma 



ma ^ — /t 3 

472T 3 fl^ -\- jm^ 27r/LL^ 



7 m 3 a 3 iii •— a 9 
%7rr dr =---: — . ri — — LA 

4&fl\ ft x + fl % 



dr 



ry«. 



(203) again. 



Another way is to calculate the variation of energy made by displacing either the 
pole m or the n% mass. The potential energy is expressed by 



\ (P 4" p) m = \ Vm + -J 2Pff/<, . 



G a 



. (205) 



where P = m/^Tr^r and ^ = 2 <r/4irr 9 the potentials of matter m//x x and <x, where r 
is the distance from m or from <r to the point where P and p are reckoned. 
The value of the second part in (205), depending upon a-, comes to 



1 H 

2 " 



fH 



m 



2 



jli-^ + yttg 4^^ . 2a 



. (206) 



and its rate of decrease with respect to a expresses the repulsion between the pole 
and the /x 2 region. This gives (203) again. 

A fourth way is by means of the quasi- electromagnetic force on fictitious interfacial 
electric current, instead of matter, the current being circular about the axis of 
symmetry AB. The formula for the attraction is 



2V curl B.E, 



r o? 



a e t * 



9 * t 



. . (207) 



if E be the radial magnetic force from m in its own medium, tensor m/Air/i^. 
Here the curl of B is represented by the interfacial discontinuity in the tangential 
induction, or 

2zm fi Y — ju 2 

4ttT 3 /^ + fa 

Also the tangential component of R is / m#/47r/A 1 r 3 , Therefore the repulsion is 



2m% fJb± 



A*3 



47rr 3 /Ltj + yu 3 47Tyti 1 r 3 



2?rr ar =. 



/*i — ^2 






47r / a 1 ^ + /t 3 



a* 



2 ., __. 



1 



a 



4^/*! fa -j- /^ 3 4;r 



(208) 



as before, equation (203). This method (207) is analogous to (204). 

§ 37* There are several other ways of representing the attraction, employing fictitious 
matter and current ; but now let us change the method, and observe how the attrac- 
tion between the magnetic pole and the iron mass is accounted for by a stress dis- 
tribution, and its space- variation. The best stress is the third, equation (188),. § 31. 
Applying this, we have simply a tension of magnitude ^/Xilii 2 = Tj in the first medium 
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and -^/x 2 R 2 3 = T 3 in the second, parallel to R x and R a respectively, each combined 
with an equal lateral pressure, so that the tensor of the stress vector is constant. 

But, so far as the attraction is concerned, we may ignore the stress in the second 
medium altogether, and consider it as the 2P N of the stress-vector in the first medium 
over the surface of the second medium. The tangential component summed has zero 
resultant ; the attraction is therefore the sum of the normal components, or 21\ cos 20 l9 
where X is the angle between R t and the normal. This is the same as 2|/x 1 (R N 2 — Rt 3 )> 
if R N and R r are the normal and tangential components of E L ; or 



for 



a 



27rr dr 1/^ 



ma 



2 H 



mz 



2^ 



A^W 3 /*j + HI \4w 3 (/*! + fa)/ 



(209) 



which on evaluation gives the required result (203). 

But this method does not give the true distribution of translational force due to the 
stresses. In the first medium there is no translational force, except on the magnet. 
Nor is there any translational force in the second ju, 2 medium. But at the interface, 
where /*, changes, there is the force — -|R 2 Vju, per unit volume, and this is represented 
by the stress-difference at the interface. It is easily seen that the tangential stress- 
difference is zero, because 

T sin 20 = /iR N R T , (210) 

and both the normal induction and the tangential magnetic force are continuous. 
The real force is, therefore, the difference of the normal components of the stress- 
vectors, and is, therefore, normal to the interface. This we could conclude from the 
expression — ^R 2 V/x. But since the resultant of the interfacial stress in the second 
medium is zero, we need not reckon it, so far as the attraction of the pole is concerned. 
The normal traction on the interface, due to both stresses, is of amount 



t iW 



/to ™ }l x 



■8*v (,1, + H y 



-\- a 



2 H ~ f-i 



(211) 



per unit area. Summed up, it gives (203) again. 

That (211) properly represents the force — * ^R 2 V/x when /x is discontinuous, we may 
also verify by supposing /x to vary continuously in a very thin layer, and then proceed 
to the limit. 

The change from an attraction to a repulsion as /x 2 changes from being greater 
to being less than /x 1? depends upon the relative importance of the tensions parallel to 
the magnetic force and the lateral pressures operative at different parts of the 
interface. In the extreme case of /x 2 = 0, we have E x tangential, with, therefore, a 
pressure everywhere. For the other extreme, R 1 is normal, and there is a pull on the 
second medium everywhere. When /x 3 is finite there is a certain circular area on the 
interface within which the translational force due to the stress in the medium 
containing the pole m is towards that medium, whilst outside it the force is the other 
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way. But when both stresses are allowed for, we see that when /x 2 >/x 1 the pull 
is towards the first medium in all parts of the interface, and that this becomes a push 
in all parts when /x x > /x 2 . 



A definite Stress only obtainable by Kinetic Consideration of the Circuital Equations 

and Storage and Flux of Energy. 

§ 38. We see that the stress considered in the last paragraph gives a rationally 
intelligible interpretation of the attraction or repulsion, The same may be said of 
other stresses than that chosen. But the use of Maxwell's stress, or any stress 
leading to a force on inductively magnetised matter as this stress does, leads us into 
great difficulties. By (198) we see that there is first a bodily force on the whole of the 
jic 2 medium, because it is magnetised, unless /x 2 =1. When summed up, the resultant 
does not give the required attraction. For, secondly, the /x x medium is also magnetised, 
unless fjb L = 1, and there is a bodily force throughout the whole of it. When this is 
summed up (not counting the force on the magnet), its resultant added on to the 
former resultant still does not make up the attraction (i.e., equivalent^ the force on 
the magnet). For, thirdly, the stress is discontinuous at the interface (though not in 
the same manner as in the last paragraph). The resultant of this stress-discontinuity, 
added on to the former resultants, makes up the required attraction. It is unneces- 
sary to give the details relating to so improbable a system of force. 

Our preference must naturally be for a more simple system, such as the previously 
considered stress. But there is really no decisive settlement possible from the theo- 
retical statical standpoint, and nothing short of actual experimental determination of 
the strains produced and their exhaustive analysis would be sufficient to determine 
the proper stress-function. But when the subject is attacked from the dynamical 
standpoint, the indeterminateness disappears. From the two circuital laws of variable 
states of electric and magnetic force in a moving medium., combined with certain dis- 
tributions of stored energy, we are led to just one stress-vector, viz. (136). It is, in 
the magnetic case, the same as (188); that is, it reduces to the latter when the 
medium, is kept at rest, so that J and Gr become J and Gr. 

It is of the simple type in case of isotropy (constant tensor), but is a rotational 
stress in general, as indeed are all the statically probable stresses that suggest 
themselves. The translational force due to it being divisible conveniently into (a) 
the electromagnetic force on electric current, (b) the ditto on the fictitious electric 
current taking the place of intrinsic magnetisation, (c) force depending upon space- 
variation of jx ; we see that the really striking part is (b). Of all the various ways of 
representing the forcive on an intrinsic magnet it is the most extreme. The magnetic 
" matter " does not enter into it, nor does the distribution of magnetisation ; it is 
where the intrinsic force h has curl that the translational force operates, usually on 
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the sides of a magnet. From actual experiments with bar magnets, needles, &c., one 
would naturally prefer to regard the polar regions as the seat of translational force. 
But the equivalent forcive 2j B has one striking recommendation (apart from the 
dynamical method of deducing it), viz., that the induction of an intrinsic magnet is 
determined by curl fa , not by h itself; and this, I have shown, is true when h is 
imagined to vary, the whole varying states of the fluxes B, D, C due to impressed 
force being determined by the curls of e and h , which are the sources of the dis- 
turbances (though not of the energy). 

The rotational peculiarity in eolotropic substances does not seem to be a very 
formidable objection. Are they not solid ? 

As regards the assumed constancy of /x, a more complete theory must, to be correct, 
reduce to one assuming constancy of /x, because, as Lord Rayleigh # has shown, the 
assumed law has a limited range of validity, and is therefore justifiable as a preparation 
for more complete views. Theoretical requirements are not identical with those of 
the practical engineer. 

But, for quite other reasons, the dynamically determined stress might be entirely 
wrong. Electric and magnetic "force" and their energies are facts. But it is the 
total of the energies in concrete cases that should be regarded as the facts, rather 
then their distribution ; for example, that, as Sir W. Thomson proved, the " mechanical 
value " of a simple closed current C is ^LC 2 , where L is the inductance of the circuit 
(coefficient of electromagnetic capacity), rather than that its distribution in space is 
given by ^HB per unit volume. Other distributions may give the same total amount 
of energy. For example, the energy of distortion of an elastic solid may be expressed 
in terms of the square of the rotation and the square of the expansion, if its boundary 
be held at rest; but this does not correctly localise the energy. If, then, we choose 
some other distribution of the energy for the same displacement and induction, we 
should find quite a different flux of energy. But I have not succeeded in making any 
other arrangement than Maxwell's work practically, or without an immediate intro- 
duction of great obscurities. Perhaps the least certain part of Maxwell's scheme, 
as modified by myself, is the estimation of magnetic energy as i|HB in intrinsic 
magnets, as well as outside them, that is, by ■JBft" ,1 B, however B may be caused. Yet, 
only in this way are thoroughly consistent results apparently obtainable when the 
electromagnetic field is considered comprehensively and dynamically. 



* 'Phil. Mag.,' Janimrv. 1887. 
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Appendix. 

Received June 27, 1891. 

Extension of the Kinetic Method of arriving at the Stresses to cases of Non-linear 
Connection between the Electric and Magnetic Forces and the Fluxes. Preserva- 
tion of Type of the Flux of Energy Formula. 

§ 39. It may be worth while to give the results to which we are led regarding the 
stress and flux of energy when the restriction of simple proportionality between 
" forces ;? and " fluxes," electric and magnetic respectively, is removed. The course to 
be followed, to obtain an interpretable form of the equation of activity, is sufficiently 
clear in the light of the experience gained in the case of proportionality. 

First assume that the two circuital laws (89) and (90), or the two in (98), hold 
good generally, without any initially stated relation between the electric force E and 
its associated fluxes C and D, or between the magnetic force H and its associated 
fluxes K and B. When written in the form most convenient for the present appli- 
cation, these laws are 



curl (H - ho) = J = C + ^ + (l) div q - BV.q), . (212) 

ot 

- curl (E - e ) = G = E + . E + (B div q - BV.q) (213) 

ot 

Now derive the equation of activity in the manner previously followed, and arrange 
it in the particular form 

e J fl + h (j -f conv Y (E — e ) (H — li ) 

= (EC + HE) + ( E ^ + H d ~) + (E J)V.q - ED div q) + (H.B V.q - HB div q), . (214) 

\ ot ot/ 

which will best facilitate interpretation. 

Although independent of the relation between E and D, &c. ? of course the dimensions 
must be suitably chosen so that this equation may really represent activity per unit 
volume in every term. 

Now, guided by the previous investigation, we can assume that (e J + h ^o) 
represents the rate of supply of energy from intrinsic sources, and also that 
V (E — e ) (H — h ), which is a flux of energy independent of q 5 is the correct form in 
general. Also, if there be no other intrinsic sources of energy than e , h 0? and no 
other fluxes of energy besides that just mentioned except the convective flux and that 
due to the stress, the equation of activity should be representable by 
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(e J + hoOo) + conv [ y ( E ~ e o) (H ~ ho) + q (U + T)] 
= (Q + tr + T) + Eq + conv %% 

= (Q + U + T)'+2av^, (215) 



where Q is the conjugate of the stress vector, F the translational force, and Q, TJ, 
and T the rate of waste and the stored energies, whatever they may be. 
Comparing with the preceding equation (214), we see that we require 

/du 3d\ /9t 3b n 

+ [E.DV.q - (ED - U) div q] + [H.BV.q - (HB - T) diy q] (216) 



Now assume that there is no waste of energy except by conduction ; then 

Q = EC + HK (217a) 

Also assume that 

3u 3d 3t 3b 

~&- B aT' dt-^Jt (2m) 

These imply that the relation, between E and D is, for the same particle of matter, 
an invariable one, and that the stored electric energy is 

.... (218) 



• »»»0»«*» 



U = 
where E is a function of D. Similarly, 

B HdB, (219) 



T 



expresses the stored magnetic energy, and H must be a definite function of B. 
On these assumptions, (216) reduces to 

2<lVg = [E.DV.q - (ED - U) diy q] + [H.BV.q - (HB - T) diy q], . . . (220) 

from which the stress-vector follows, namely, 

F N = [E.DN - N (ED - U)] + [H.BN - N (HB ~~ T)]. (221) 

Or, 

P N = (VDVEN + NU) + (VBVHN + NT). (222) 

Thus, in case of isotropy, the stress is a tension U parallel to E combined with 
a lateral pressure (ED — U) ; and a tension T parallel to H combined with a lateral 
pressure (HB — T). 

The corresponding translational force is 

F = E diy D + DV.E - V (ED - U) 

■f HdivB + BV.H- V(HB-T), (223) 

which, it is unnecessary to put in terms of the currents. 
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Exchange E and B, and H and B, in (221) or (222) to obtain the conjugate vector 
Q N ; from which we obtain the flux of energy due to the stress, 

- qQ, q = D.Eq - q (ED - U) + B.Hq - q (HB - T) 

= YlYBq 4- YHYBq + q (U + T), (224) 

or 

- ^ = YeH + YEh + q (U + T), (225) 

where e and h are the motional electric and magnetic forces, of the same form as before 
(88) and (91); so that the complete form of the equation of activity, showing the 
fluxes of energy and their convergence, is 

e^ +hoG +conv[V(E-e )(^^ (226) 

where F has the ahove meaning. 

There is thus a remarkable preservation of form as compared with the corresponding 
formulae when there is proportionality between force and flux. For we produce 
harmony by means of a Poynting flux of identical expression and a flux due to the 
stress, which is also of identical expression, although U and T now have a more general 



meaning, of course. 



•Jii. 



Example of the above, and Remarks on Intrinsic Magnetisation when there is 

Hysteresis, 

§ 40. In the stress-vector itself (for either the electric or the magnetic stress) the 
relative magnitude of the tension and the lateral pressure varies unless the curve 

# As the investigation in this Appendix has some pretensions to generality, we should try to settle 
the amount it is fairly entitled to. ISTo objection is likely to be raised to the use of the circuital equations 
(212), (213), with the restriction of strict proportionality between E and H and the fluxes D and B, or 
C and K entirely removed ; nor to the estimation of J and Gr as the " true " currents ; nor to the use of 
the same form of flux of electromagnetic energy when the medium is stationary, For these things are 
obviously suggested by the preceding investigations, and their justification is in their being found to 
continue to work, which is the case. But the use in the text of language appropriate to linear functions, 
which arose from the notation, &c, being the same as before, is unjustifiable. We may, however, remove 
this misuse of language, and make the equation (226), showing the flux of energy, rest entirely upon 
the two circuital equations, In fact, if we substitute in (226) the relations (217a), (2176), it becomes 
merely a particular way of writing (214). 

It is, therefore, to (217a), (2176) that we should look for limitations. As regards (217a), there does 
not seem to be any limitation necessary. That is, there is no kind of relation imposed between E and C> 
and H and IL This seems to arise merely from Q meaning energy wasted for good, and havings no 
further entry into the system. But as regards (2176), the case is different* For it seems necessary, in 
order to exclude terms corresponding to E(3c/3£)E and H(3/*/3j^)H in the linear theory, when there is 
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connecting the force and the induction be a straight line. Thus, if the curve be 
of the type shown in the first figure, the shaded area will represent the stored 
energy and the tension, and the remainder of the rectangle will represent the lateral 
pressure. They are equal when H is small ; later on the pressure preponderates, and 
more and more so the bigger H becomes, 




But if the curve be of the type shown in the second figure, then, after initial 
equality the tension preponderates ; though, later on, when H is very big the 
pressure preponderates. 

To obtain an idea of the effect, take the concrete example of an infinitely long rod, 
uniformly axially inductized by a steady current in an overlapping solenoid, and 
consider the forcive on the rod. Here both H and B are axial or longitudinal ; and 
so, by equation (223), the translational force would be a normal force on the surface of 
the rod, acting outwards, of amount 

— 1) — -^KqBq 



per unit area ; this being the excess of the lateral pressure in the rod over -| H B , the 
lateral pressure just outside it. 

In case of proportionality of force to flux, the first pressure is ^HB, and if there is 
no intrinsic magnetisation H and H are equal. The outward force is therefore 
positive for paramagnetic, and negative for diamagnetic substances, and the result 
would be lateral expansion or contraction, since the infinite length would prevent 
elongation. 

rotation, that E and D should be parallel, and likewise H and B. At any rate, if such terms be allowed, 
some modification may be required in the subsequent reckoning of the mechanical force. In other respects, 
it is merely implied by (217b) that E and D are definitely connected, likewise H and B, so that there is 
no waste of energy other than that expressed by Q. 
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But if the curve In the rod be of the type of the first figure, and the straight 
line ac be the air curve to correspond, it is the area abc that now represents the 
outward force per unit area when the magnetic force has the value ad. If the 
straight line can cross the curve ab, we see that by sufficiently increasing H we can 
make the external air pressure preponderate, so that the rod, after initially expand- 
ing, would end. by contracting. 

If the rod be a ring of large diameter compared with its thickness, the forcive 
would be approximately the same> viz., an outward surface force equal to the 
difference of the lateral pressures in the rod and air. The result would then be 
elongation, with final retraction when the external pressure came to exceed the 
internal. 

Bid well found a phenomenon of this kind in iron, but it does not seem possible 
that the above supposititious case is capable of explaining it, though of course the 
true explanation may be in some respects of a similar nature. But the circumstances 
are not the same as those supposed. The assumption of a definite connexion between 
H and B, and elastic storage of the energy T, Is very inadequate to represent the facts 
of magnetisation of iron, save within a small range. 

Magneticians usually plot the curve connecting H — h and B, not between H and B, 
or which would be the same, between II — h and B — B Q3 where B is the intrinsic 
magnetisation. Now when an iron ring is subjected to a given gaussage (or magneto 
motive force), going through a sequence of values, there is no definite curve connecting 
H — h Q and B, on account of the intrinsic magnetisation. But, with proper allowance 
for A , it might be that the resulting curve connecting H and B in a given specimen, 
would be approximately definite, at any rate, far more so than that connecting H — h 
and B. Granting perfect definiteness, however, there is still insufficient information 
to make a theory. The energy put into Iron is not wholly stored ; that is, in 
increasing the coil current we increase B as well as B, and in doing so dissipate 
energy ; and although we know, by Ewing's experiments, the amount of waste in 
C3 7 clical changes, it is not so clear what the rate of waste is at a given moment. 
There is also the further peculiarity that the energy of the intrinsic magnetisation at 
a given moment, though apparently locked up, and really locked up temporarily, 
however loosely it may be secured, is not wholly Irrecoverable, but comes into play 
again when H is reversed. Now it may be that the energy of the intrinsic magnetisa- 
tion plays, in relation to the stress, an entirely different part from that of the elastic 
magnetisation. It is easy to make up formulae to express special phenomena, but 
very difficult to make a comprehensive theory. 

But in any case, apart from the obscurities connected with iron, it is desirable to be 
apologetic in making any application of Maxwell's stresses or similar ones to 
practice when the actual strains produced are In question, bearing in mind the 
difficulty of interpreting and harmonising with Maxwell's theory the results of 
Kerr, Quincke, and others. 



